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ZUSAMHENFASSUNG 


In  diesem  Bericht  wird  in  Experiment  und  Theorie  die  Art 
der  Ausbreitung  von  VLF-\7ellen  im  festen  Medium  fur  das 
Nahfeld  und  die  anschlieBende  lib ergangs zone  zum  Fernfeld 
beschrieben. 

Die  maBgebenden  Parameter  &  (Leitfahigkeit)  und  £.  (Dielektri- 
zitatskonstante)  werden  durch  unabhangige  MeBverfahren 
bestimmt.  Die  notwendigen,  vielfach  selbst  konstruierten 
Gerate  werden  beschrieben.  Durch  eine  Vielzahl  von  Messungen 
in  geologisch  verschiedenen  Gebieten  wird  die  Theorie  des 
strahlenden  magnetischen  Dipols  uberpriift  und  modifiziert. 


ABSTRACT 


In  the  present  report,  the  propagation  of  VLF  waves  in  solid 
media  is  described  by  theory  and  experiment  for  the  near 
field  and  for  the  subsequent  transition  zone  to  the  far  field. 

The  decisive  parameters  6^ (conductivity)  and  £.  (dielectric 
constant)  are  determined  by  independent  methods  of  measure¬ 
ment.  The  required  instruments,  largely  designed  by  our 
team,  are  described.  The  theory  of  a  magnetic  transmitting 
dipole  is  checked  and  modified  by  a  number  of  measurements 
in  regions  of  different  geological  conditions. 
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INTRODUCTION 

The  propagation  of  electromagnetic  ?mves  in  air,  i.e., 
in  the  earth  -  ionosphere  cavity  as  well  as  outside  the  ionos¬ 
phere  has  "been  studied  sufficiently  in  the  region  from  long 
waves  to  ultrashort  waves.  Also  the  studies  on  the  propagation 
mechanism  of  very  long  (VLF)  to  extremely  long  (ELF)  waves 
in  the  atmosphere  have  almost  been  ooncluded  in  the  past  few  years. 
The  behavior  of  VLF  waves  passing  through  homogeneous  or  in¬ 
homogeneous  media  was  studied  only  as  far  as  these  problems 
were  of  geophysical  interest.  The  studied  frequencies  mainly 
lay  in  the  long  wave  region  or  short  wave  region,  i.e.,  at 
frequencies  above  50  keps. 

The  VLF  project  Dr.  W.  Bitterlich  which  was  founded  in 
autumn  I960  under  the  sponsorship  of  the  US  Government,  had 
the  aim  of  studying  the  propagation  mechanism  of  long  electro¬ 
magnetic  waves  (<  20  keps)  passing  through  the  interior  of 
the  earth.  AL  the  same  time  the  resulting  possibilities  of 
underground  communication  had  to  be  studied  as  well  as  the 
possibilities  of  their  practical  application  in  geophysics, 
geology  and  mining. 

The  present  report  gives  a  review  of  all  VLF  studies  made 
between  December  i960  and  December  1965*  The  principal  directions 
of  this  work  were  the  following: 

1.  Investigation  of  the  propagation  of  VLF  waves  in  rock  by 
means  of  our  own  transmitters;  preferred  frequencies:  5  and 
10  keps.  The  field  strengths  of  measuring  distances  given 

at  various  points  were  determined  and  evaluated  as  functions 
of  the  transmitter  -  receiver  distance  and  the  angle  of 
direction  of  the  axes  of  transmitting  antennas  and  receiving 
antennas . 

2.  Studies  of  the  penetration  of  electromagnetic  waves  of 

a  -  w 

remote  VLF  transmitters  into  the  ground. 

5.  The  directivity  pattern  of  a  magnetic  dipole  was  determined 
a3  a  function  of  conductivity  of  the  surrounding  medium  of 


the  transmitter  -  receiver  distance  and  of  the  frequency. 

4.  Determination  of  the  electrical  conductivity  of  natural 
rock  with  various  direct  current  and  alternating  current 
methods  by  direct  contact,  and  studies  on  the  frequency  de¬ 
pendence  of  the  dielectric  constant  of  rock  samples. 

.  5«  Detailed  theoretical  studies  of  problems  regarding  pos.  1-4. 

At  the  end  of  the  present  report,  the  titles  of  the  pre¬ 
vious  16  Scientific  Reports  are  given  as  well  as  the  tables 
of  contents  of  the  four  Technical  Annual  Summary  Reports. 

During  the  entire  period  of  project  work,  Dr.  0.  Grobner, 
Dr.  0  Worz,  VT.  Gradl,  G.  Tinhofer  and  N.  Nessler  were  my 
closest  collaborators.  The  above  gentlemen  contributed  de¬ 
cisively  to  the  results  and  I  herewith  thank  them  cordially 
for  their  collaboration. 

The  measurements  in  the  mines  were  made  possible  owing  to 
the  generious  support  by  Bergrat  Dipl.  Ing.  P.  Kettner, 

Director  of  the  Montanwerke  Brixlegg  who  permitted  us  the 
establishment  of  a  mine  laboratory,  and  the  friendly  assistance 
of  the  Berghauptraanns chaf t  of  the  province  of  Tyrol.  Let  me 
now  express  my  sincere  gratitude  to  the  concerned  gentlemen 
and  above  all  to  Univ.  Prof,  Dr.  R.  Steinmaurer,  Head  of  the 
Physics  Institute  of  the  University  of  Innsbruck. 


December  1965 


Dr.  Wolfram  Bitterlich 
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1.  Apparatus 

1.1.  General  Information 


At  the  beginning  of  our  studies  in  1961,  a  considerable 
number  of  instruments  that  are  now  being  produced  in  series 
were  not  on  the  market.  Especially  battery-fed  power  amplifiers 
of  low  weight  and  high  efficiency  as  well  as  mains-indepondent 
selective  level  meters  for  the  VLF  region  were  not  available. 

We  therefore  had  to  rely  on  our  own  developments  and  special 
designs.  The  development  and  testing  of  suitable,  i.e.,  small 
antennas  required  great  efforts  in  time  and  work. 

Instruments  used  for  measurements  below  ground  are  ex¬ 
posed  to  much  harder  conditions  than  those  used  in  the  laborato¬ 
ry.  Transportation  to  the  sites  of  measurement  is  troublesome, 
as  it  is  only  in  the  main  galleries  that  instruments  can  be 
transported  with  a  mine  railway,  and  even  there  they  are  exposed 
to  hard  shocks  and  rough  handling.  The  actual  region  of  measure¬ 
ment  should  be  far  away  from  tracks,  tubes  and  cables,  it  is 
therefore  only  accessible  on  foot  through  rather  difficult 
terrain.  Occasionally,  remote  sites  are  only  reached  by  means 
of  ladders,  manshafts  or  by  crawling.  The  extremely  high  water 
content  in  the  atmosphere  (100^6,  numerous  galleries  are 
dripping  wet  with  water  running  over  the  bottom  of  the  gallery) 
makes  special  demands  to  insulation  and  antioorrosion. 

Prom  these  requirements  result  the  partly  contradictory 
conditions  for  developing  and  constructing  special  instruments 
for  our  experiments.  The  devices  must  be  strong,  shockproof, 
protected  against  dripping  water  and  provided  with  a  waterproof 
base.  Yet  they  must  be  light  and  sufficiently  small  or 
collapsible  as  to  pass  through  narrow  passages  in  the  form  of 
single  parts. 
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1.2.  Transmitter 


In  order  to  make  the  efficiency  as  high  as  possible,  we 
used  transistorized  push-pull  amplifiers,  except  for  a  prelimi¬ 
nary  LF  experiment.  Fig.  1.1  shows  the  device  used  for  our 
final  measurements.  A  transistorized  RC  generator  serves  as 
an  oscillator.  The  driver  is  designed  as  usual.  The  output-stage 
in  C-operation  works  on  an  autotransformator ,  thus  keeping  the 
transfer  losses  low  and  yet  offering  a  good  possibility  of 
antenna  adaptation.  The  attained  compensating  power  was  50  w, 
the  cut-off  frequency  was  10  kcps.  The  antenna  current  is 
indicated  on  a  built-in  instrument  via  a  control  transformer. 

For  tuning  the  antenna,  a  multistage  capacitor  is  built  in 
which  consists  of  mica  capacitors  with  a  maximum  dielectric 
strength  of  2000  volts.  In  order  to  protect  the  power  transistors 
against  voltage  peaks  that  may  occur  at  the  moment  when  the 
antenna  is  switched  over,  the  collectors  are  short-circuited 
with  powerful  silicon  diodes  as  limiters. 

Because  of  the  weight  restriction  of  the  individual 
parts,  the  available  battery  capacitance  is  limited  to  about  75  ah. 
\7ith  an  output  of  approximately  50  w  and  an  efficiency  of  about 
50$  we  have  a  working  time  of  7  hours  which  is  quite  sufficient 
for  long-time  experiments  (description  see  \j>j  ). 

For  attaining  a  higher  transmitter  power,  we  have  set 
up  a  1  kw  -  IF  amplifier  (Bryan  Savage)  type  KM2Z  in  our 
laboratory  in  the  mine  cf  St.  Gertraudi.  The  amplifier  has  a 
frequency  ranging  from  20  cps  to  10  kcps  at  full  output.  The 
matching  impedance  can  be  chosen  between  0.6  and  40  £2. 

In  order  to  extend  the  frequency  range  toward  higher 
values  we  used  a  stereoamplifier  (McIntosh)  type  MC275  for 
our  recent  measurements.  It  has  an  output  of  150  w  up  to  40  kcps 
and  75  w  up  to  100  kcps,  with  change-over  output  impedances  of 
4  to  600  Q. 
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1.3.  Transmitting  antennas 

Based  upon  practical  and  theoretical  considerations, 
especially  with  respect  to  power  adaptation  and  dimensions, 
we  have  used  magnetic  dipoles  as  transmitting  antennas  right 
from  the  beginning  (  ! ^chapter  1.2  and  1.3,  f  3 1  chapter  2.3, 
fed).  At  first,  rod  antennas  with  iron  cores  made'  of  trans¬ 
former  sheet  metal  with  center-tapped  windings,  were  connec¬ 
ted  directly  to  the  power  transistors  of  the  push-pull  stage 
without  output  transformer.  The  magnetic  moment  of  this 
array,  however,  was  too  small  because  of  the  small  turn-area. 
Tfe  therefore  used  air  coils  with  possibly  large  diameter, s 
instead  of  the  above  coils.  As  described  in  Fl}  ,  they  were 
chosen  in  accordance  with  the  given  power  of  amplifier  and 
the  possible  weights  and  dimensions. 

Data  of  the  transmitting  antennas  that  were  specially 
suited  for  propagation  measurements: 

SA  VI  [2]  square  frame,  area  1.4  m? 

2 

60  turns  of  copper-stranded  wire,  2.5  mm  cross 
section,  15  mH,  2.1  il. 

In  this  antenna,  the  copper-stranded  wire  is  welded 
as  a  bunch  in  a  flexible  plastic  hose.  Thus,  the 
winding  can  be  removed  from  the  frame  and  folded  up, 
the  frame  itself  being  collapsible.  Transportation 
of  the  antenna  is  thus  very  easy. The  majority  of 
measurements  with  this  antenna  were  made  by  means 
of  a  battery-fed  transmitter,  above  all  those  in 
the  German  and  Dutch  mines  (cf.  chapter  2). 

SA  VII 1 3]  Rectangular  frame,  area  17  m^. 

2 

2  x  15  turns,  1.5  mm  copper-stranded  wire,  12  mH, 
13.5  This  antenna  is  suspended  in  the  "Bunte 
Kluft"  in  the  St.  Gertraudi  mine  and  can  be  rotated 
about  its  vertical  axis.  It  has  been  used  for 
numerous  measurements  with  the  1  kw  transmitter. 
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SA  VIII  [ 3~1  An  approximately  square  loop  having  an  area  of 

p  2 

1600  m‘,  3  turns,  2.5  mm  copper,  20  mH,  3  Q, 

In  order  to  attain  these  dimensions,  the  antenna 
had  to  be  firmly  established  at  a  suitable  site 
at  St.  Gertraudi  in  two  shafts  with  a  vertical 
distance  of  approximately  40  m  and  in  two  oonnecting 
galleries  having  a  vertical  interval  of  4°  m. 

SA  IX  As  described  in  £4]  >  the  SA  VIII  had  been  re¬ 

constructed.  The  number  of  turns  was  increased 

2 

to  10  and  the  cross  section  to  135  mm  (aluminum). 
Thus  it  was  possible  to  get  antenna  currents  of 
up  to  40  a. 

The  transmitting  antennas  were  tuned  to  series  resonance 
so  that  the  power  can  be  reached  by  means  of  the 

tuning  capacitor,  thus  only  the  active  power  (which  almost 
completely  was  consumed  in  the  wire  in  the  form  of  ohmic  losses) 
being  attained  by  the  transmitter.  Because  of  the  necessary 
alternating  voltage  strength  and  bhe  high  reactive  currents, 
mica  capacitors  were  used  for  lower  powers  and  oil  capacitors 
for  the  1  kw  amplifier. 


In  order  to  keep  the  dimensions  and  weights  as  low  as 
possible,  only  transistorized  receivers  were  developed  [l  - 
Figs.  1.2  -  1.5  show  the  circuit  of  a  highly  suited  in¬ 
strument.  The  input  stage  is  equipped  w ith  two  low-noise 
transistors  in  emitter  circuit  and  a  silicon  transistor  in 
collector  circuit,  serving  as  impedance  transformer  for  the 
stage  attenuator.  The  input  voltage  on  the  first  transistor 
is  limited  to  a  maximum  of  approximately  100  mv  by  means  of 
the  OA  71  diodes  in  order  to  avoid  overloading  of  the  input 
stage.  The  range  of  measurement  is  switched  over  with  S  1 
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such  that  attenuation  before  the  input  stage  is  4O  <13  and 
behind  it  10  dB.  The  input  resistance  is  kept  constant  by  the 
switch  plane  Sib.  The  second  attenuator  is  connected  with  a 
selective  amplifier.  A  bridged  T-type  section  has  a  frequency- 
dependent  negative  feed-back  with  three  stages.  Tuning 
is  performed  in  stages  1  :  10  :  100  by  switching  over  the 
capacitors,  and  continuously  1  :  10  by  a  10- turn  helical 
potentiometer  with  two  sections.  The  band  width  can  be  ad¬ 
justed  by  changing  a  feedback  which  is  independent  of  the 
frequency.  At  sufficient  stability,  a  relative  bandwidth  of 
=  5  -  10 °/o  was  attained.  The  selective  amplifier  is  followed 
by  an  indicating  amplifier  feeding  an  instrument  which  is 
calibrated  directly  in  pv  or  mv.  For  controlling  the  voltage 
calibration,  a  calibration  oscillator  is  built  in  whose 
calibration  voltage  stabilized  with  a  diode  circuit  is 
sufficiently  constant  at  temperatures  between  +5°C  and  +  35°C. 
A  superheterodyne  supplement  permits  the  acoustic  control  of 
nonmodulated  or  keyed  transmitters.  The  operating  voltage 
of  most  stages  is  stabilized  electronically. 

Technical  data:  Measuring  range  1  pv  to  100  mv  full  deflection 

accuracy  of  measurement  +  5 c!° 
noise  voltage  of  input  0.1 
input  resistance  1.9  kQ 
frequency  range  20  cps  to  30  kcps 
relative  band  width  5 
detailed  description  in  |  3_|  • 

Nuvistorized  field  strength  meter: 

For  further  improvement  of  the  sensitivity  and  the 
signal-to-noise  ratio,  nuvistors  were  substituted  for  the 
transistors.  Thus,  a  very  high  input  resistance  was  attained 
which  on  the  one  hand  kept  the  high  quality  of  a  ferrite  an¬ 
tenna  tuned  to  parallel  resonance  as  well  as  its  high  se¬ 
lectivity  and  induced  voltage,  and  on  the  other  hand  permits 
the  useful  application  of  high  quality  resonant  circuits  in 
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the  subsequent  circuits  too. 

An  input  stage  according  to  Pig.  1.6,  an  improved  version 
of  [12],  fixed  to  the  antenna.  It  consists  of  a  cascode 
connection  with  a  triode  and  a  pentode  followed  by  another 
triode  as  impedance  transformer.  The  antenna  coil  is  symmetri¬ 
cal,  its  center  tap  being  grounded.  Coupling  with  the  non- 
symmetrical  cascode  input  is  made  with  a  transducer.  Antenna 
and  transducer  are  jointly  tuned  to  resonance.  The  device 
contains  a  calibration  oscillator  for  adjusting  the  amplifi¬ 
cation  by  means  of  the  potentiometer  R  1.  Low-ohmic  output 
voltage  is  collected  on  Bu2  and  further  amplified  and  measured 
by  separate  instruments.  Current  supply  is  provided  either 
jointly  with  our  own  field  strength  meter  discussed  in  [12]], 
or  separately  with  a  6  v  lead  battery  and  a  120  v  nickel- 
cadmium  battery.  Following  such  an  input  stage,  the  selective 
measuring  amplifier  described  ir.  detail  in  £l2l  was  built.  It 
is  also  nuvistorized  using  an  LC  series  resonant  circuit 
as  a  means  of  selection  which  is  connected  between  the  output 
cathode  follower  of  the  input  stage  and  the  input  of  the  meter 
amplifier  such  that  the  considerable  voltage  increase  due  to 
the  high  quality  is  fully  utilized.  Technical  data  of  the 
input  stage  of  Fig.a6: 

Maximum  input  voltage  on  Bu  1:  200  mv 
Amplification  factor:  50 

noise  on  the  short-circuited  antenna  input:  0.2  V-Vef.£«> 
frequency  range:  2.4  -  26  keps, 

band  width  of  the  antenna  FA  7  with  input  stage  alone 
at  10  keps:  75  cps,  at  5  keps:  10  cps. 

Together  with  the  transmission  level  meter  Marcony  TF  2350  a 
maximum  sensitivity  range  of  0.6  pv  full  deflection  is  reached. 
Besides  our  own  meters,  commercial  devices  were  used  in  the 
course  of  our  studies.  We  used  the  fallowing  instruments: 
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Wave  Analyzer  TP  2330,  Marconi  Instruments  Ltd,  Great  Britain 
Frequency  range  20  cps  -  $0  kcps. 

Measuring  range  30  -  300  v  full  deflection, 

Measuring  accuracy  5$> 

Input  resistance  100  kQ/lMQ 
Band  width:  3.5  cps, 

Battery-fed. 

Y/ave  Analyzer  KLA  48,  Wandel  und  Goltermann,  Germany 
Frequency:  30  cps  -  50  kcps, 

Measuring  range:  0.25  mv  -  8v  full  deflection 
Band  width  changeable,  5*5  and  15  cps,  re¬ 
spectively. 

Wave  Analyzer  SPM  2,  Wandel  void  Goltermann,  Germany 
Frequency  region  4  -  600  kcps , 

Band  width:  200  cps, 

Measuring  range  2  mv  »  7  v  full  deflection. 
Battery-fed. 

Oscilloscope  421,  Tektronix,  USA. 

1.5  Receiving  antennas 


It  has  already  been  shown  in  |  4 1  and  fl4^[  that  preferably 
magnetic  receiving  antennas  should  be  used  for  our  present 
problems,  mainly  because  of  their  easy  transportation.  For 
the  induced  voltage  Ue  on  a  coil  tuned  to  resonance,  the 
following  relation  is  valid: 


Ue=w.n.F.B.v.effQ.k 


U) 

n 

F 

^eff 

k 


angular  frequency 
number  of  turns 

turn  area  (in  air  coils)  or  core  cross  section  (ferrite 
ruds ) 

effective  permeability  for  the  used  core  shape 
form  factor  which  depends  on  the  coil  array,  lies  bet¬ 
ween  approximately  0.5  and  1. 


O 
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In  air  coils,  =  1  and.  the  induced  voltage  can  only  be 

increased  by  increasing  the  turn  area  nP.  Apart  from  the  de¬ 
sired  restriction  in  size,  the  co -.1  losses  thus  increase  con¬ 
siderably  so  that  the  reduction  of  Q  obstructs  the  de¬ 
sired  effect.  By  means  of  red  cores  of  highly  permeable 
ferrites  which  are  well  suited  for  VLF,  a  p  ^  of  about  50 
to  100  can  be  attained.  A  high  permeability  involves  high 
inductance  and  thus  also  a  large  L/C  ratio.  Together  with 
tile  small  diameter  of  winding  which  yields  short  wire  length 
and  small  coil  losses,  a  high  Q  (usually  above  100)  is  ob¬ 
tained  which  can  only  be  used  if  the  input  resistance  of  the 
connected  amplifier  input  circuit  is  correspondingly  high 
as  for  example  in  the  field  strength  meter  E4  described  in  [l2j  . 

Of  all  tested  antennas,  the  most  typical  and  preferred 
types  are  described  hereinafter: 

EA  I  :{' Ijferrite  rod,  length  47.5  cm,  diameter  1.5  cm,  in¬ 
ductance  0.8  H,  ohmic  resistance  155 
EA  III: [ lj cylinder  coil  without  iron  core,  diameter  45  cm, 

1000  turns. 

Both  antennas  are  tuned  to  parallel  resonance,  the  receiver  fl| 
(10  k&  input  resistance)  is  always  connected  with  a  tap.. 

At  first,  the  antenna  EA  III  was  used  for  numerous  measure¬ 
ments  . 

EA  IV:  [2]  five  ferrite  rods  of  2  cm  diameter  and  80  cm 

length  each,  arranged  as  the  edges  of  a  regular  penta¬ 
gonal  prism.  5  x  1600  turns,  5  x  50  Q,  always  in  the 
rod  center;  total  inductance  has  a  maximum  of  8  H. 

EA  V:  [2] a  ferrite  rod  of  2  cm  diameter  and  65  cm  length, 

5000  turns  in  the  middle,  155  2  H. 

EA  VI : [2J a  laminated  iron  core  made  of  metal  sheet  strips 

(0.55  nim  thick  and  50  mm  wide),  total  iron  cross  section 
2 

8  cm  ,  length  100  cm,  material  Hyperm  7 66  by  Krupp-Widia, 

Essen,  with  a  very  high  ring  core  permeability  cf  u,  = 

tor 

=  80000. 
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The  improved  transistor  receivers  had  a  smaller  input 
resistance  (E  3:  1.9  kQ) .  The  antennas  were  operated  in 
series  resonance.  At  the  receiver  input  E  2  [2]  they  yielded 
a  voltage  of  1  yv  at  a  field  strength  of 

_ f  =  3  kcps _ f  =  10  kcps _ 

EA  IV  2.32  •  I0~]l  O.73  •  10"^>7b.m”2 

ea  v  6.8  •  10":^  2.42  .  10"r^  " 

EA  VI  3.25  •  10  2.63  •  10"  d  « 

The  input  voltage  of  1  yv  at  the  receiver  E  3  [24] 

was  reached  at: 

_ f  =  3  kcps _ f  =  10  kcps 

EA  IV  0.26  •  10~\2  0.15  •  lO'^Tb.nf 2 

EA  VI  1.3  .  10-1^  0.55  •  10  " 

FA  V:  Ferrite  rod  of  2  cm  diameter  and  48  cm  length:  in  the 
middle  1200  turns,  450  mH. 

This  antenna  is  installed  as  a  direction-finder  antenna 
in  Cardanic  suspension  so  that  the  value  as  well  as  the 
direction  of  the  field  strength  can  he  measured  accurately. 

Fig.  1.9  shows  the  principle  of  the  array  that  can  he  rotated 
in  all  three  space  axes.  In  order  to  avoid  a  field  distortion 
by  metal  parts,  the  construction  was  entirely  made  of  plastics. 
The  antenna  is  also  tuned  to  parallel  resonance  and  operates 
on  an  impedance  transformer  which  although  having  a  voltage 
amplification  <  1,  permits  the  high  quality  of  the  antenna 
circuit  to  he  utilized.  The  impedance  transformer  (Fig.  1.10) 
consists  of  two  collector  stages  in  cascade,  equipped  with 
HF  transistors.  A  considerable  reduction  of  pickup  from 
electrical  stray  fields  was  attained  by  tapping  and  grounding 
the  coil  center  of  the  antenna  and  by  connecting  the  impedance 
transformer  symmetrically  via  a  transformer  with  symmetrical 
input  and  ronsymrnetrical  output.  The  transformer  transforms 
upward  at  1  :  1.8  and  is  tuned  to  resonance  together  with  the 
antenna.  At  10 kcps,  this  improved  antenna  yields  a  voltage  of 
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1  pv  at  7*75  '  10*"^  ,Tb*rn”^  at  the  output  of  the  impedance 
transformer.  The  calibration  curves  for  different  frequencies 
and  field  strengths  are  shown  in  Pig.  1.8. 

FA  7s  Ferrite  rod  of  2  cm  diameter  and  98  cm  length.  1200  turns 
in  the  middle. 

This  antenna  was  originally  developed  for  the  E  4  re¬ 
ceiver,  for  details  see  [l2 J  and  [14I .  Together  with  the 
improved  input  stage  as  shown  in  Fig.  1.6,  the  FA  7  was 
changed  in  the  same  way  as  FA  5  with  respect  to  symmetrical 
winding  and  transformer  connection.  At  the  output  of  the 

first  stage  (Bu  2,  Fig.  1.6),  FA  7  yields  a  voltage  of  1  mv 
-15  -2 

at  6.6  •  10  J  Vfb*m  at  10  kcps.  Fig.  17  shows  the  calibration 
curves  as  dependent  on  the  frequency.  Here,  the  increase  in 
LC  ratio  at  higher  frequencies  becomes  especially  evident. 

The  discontinuity  at  5  kcps  corresponds  to  the  impedance 
switch-over  with  S  1  of  Fig.  1.6. 

For  the  antennas  EA  IV  and  EA  V  as  well  as  FA  5  and  FA  7> 
the  material  1100  N22  of  Siemens  und  Haiske,  Munich  was  used. 

1.6  Calibration  of  field  strength  meters 

As  the  voltage  induced  in  the  antenna  depends  on  the 

frequency,  calibration  curves  of  the  used  antennas  together 

with  the  receivers  were  made  as  a  function  of  frequency. 

For  this  purpose,  a  calibration  field  was  produced  according 

to  Biot-Savart  ’  s  lav/,  A  circular  calibration  coil  with  an 
2 

area  of  1290  cm  and  10  turns  was  fed  by  an  RC  generator,  the 
current  being  measured  with  a  clip-on  a-c  current  probe. 

A  series  resistor  of  150  Q  betv/een  the  generator  and  the 
coil  inductance  of  75  pH  made  the  coil  current  sufficiently 
independent  of  the  frequency  up  to  about  20  kcps,  up  to 
which  frequency  the  maximum  impedance  of  the  coil  is  no  more 
than  l/l5  of  the  chosen  value  of  the  barrier  resistance.  The 


f  [kHi] 


Parameter  I  Hi  [  Wh/rn1'] 


Fig..  1,6 
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receiving  antenna  to  be  calibrated  was  set  up  at  a  sufficient 
distance  which  is  well  known,  coaxially  to  the  calibration  coil. 
Under  the  condition  that  the  field  distortion  due  to  the 
introduction  of  a  receiving  antenna  is  negligible,  the 
following  relation  holds  for  the  field  strength  at  the  site 
of  the  receiving  antenna  : 

2 

H  =  2npo*I*n* 

P 

H  =  field  strength  in  \!b/ m  : 
n  =  number  of  turns  of  the  calibration  coil 

r  =  distance  between  the  plane  of  the  calibration  coil  and  the 
end  of  the  receiving  antenna  facing  the  calibration  coil,  m, 
I  =  coil  current,  a, 
a  =  radius  of  the  calibration  coil. 

In  general,  the  radius  a  is  small  as  compared  to  the 
distance  r,  therefore  the  expression  is  simplified  as  follows: 

2u  •  I  •  n  •  F 


For  practical  reasons,  the  distances  used  for  the  calibration 
varied  between  5  and  15  m,  the  coil  current  reached  about  6  ma. 
Used  instruments: 

Decade  PC  generator  TG  66,  Levell,  England 

Clip-on  a-c  current  probe  456  A,  Hewlett-Packard  Company,  USA. 
Calibration  v/as  controlled  by  small  receiving  calibration 
coils  described  in  j^2]  ,  to  guarantee  that  the  field  change  of 
the  calibration  field  caused  by  introducing  the  receiving  an¬ 
tennas  was  definitely  kept  below  the  admissible  limits. 

Figs.  1.7  and  1.8  show  examples  of  the  calibration  curves 
obtained  in  the  above  manner. 


f  2  2s 

(r  +  a  ) 


2-1 


2.  Experimental  studies  in  the  mine 

2.1  Introduction 

Measurements  below  ground  were  made  for  the  purpose 
of  determining  the  manner  of  VLF-wave  propagation  in  a  solid 
medium  and  the  effect  and  quantity 

of  the  parameters  £. ,  p,  S'  (dielectric  constant,  per¬ 
meability,  conductivity)  in  connection  with  it,  as  well  as 
the  relation  between 

the  transmitter  power  (value  of  the  magnetic  moment  of 
the  transmitting  coil  m)  and 

the  transmitter  position  (angle  -J-  between  the  trans¬ 
mitting  dipole  and  the  radius  vector  with  respect  to  the 
point  of  measurement), 

distance  between  transmitter  and  receiver  and  th^ 
angularity  of  the  receiving  antenna  (angle  y  ). 

The  permeability  p  with  good  approximation  was  found 
to  be  un.  ity  for  all  types  of  rock.  Measurements  showed  this 
to  be  valid  with  sufficient  accuracy  also  for  hematitic  and 
sideritic  ores.  The  dielectric  constant  t  was  measured  on 
cylindrical  laboratory  samples.  The  results  are  summarized 
in  a  special  chapter  of  the  present  report.  On  measuring 
the  conductivity  S'  it  was  soon  found  out  that  the  conditions 
of  measurement  influence  the  results  in  a  rather  complex 
way.  £  was  therefore  determined  by  two  completely  separate 
methods:  by  direct  measurement  on  solid  rock  by  the  Wenner 
method  (see  section  5)  and  by  measuring  the  propagation  of 
VLP  waves. 

The  propagation  measurements  were  made  in  different 
ways,  in  accordance  with  the  numerous  problems  that  arose. 

It  has  already  been  said  above  that  the  transmitter 
(voltage  source,  generator  of  accurate  frequency,  amplifier 
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and  transmitting  antenna)  were  set  up  in  a  suitable  cavern 
of  the  mine,  the  transmitting  antenna  being  aligned  under 
a  defined  angle  by  means  of  a  geological  compass.  In  most 
cases,  the  axis  of  the  transmitting  frame  and  thus  also 
of  the  magnetic  dipole  has  a  horizontal  position.  The 
frame  (in  rotating  antennas)  is  moved  about  a  perpendicular 
axis  (angle  ,  see  Pig.  2.l). 

The  measurement  can  thus  be  made  in  the  same  manner 
as  a  profile  measurement,  i.e.,  the  field  strength  is 
measured  by  means  of  the  receiver  as  rectilinear ly  as  possible, 
at  suitable  distances  from  the  transmitter.  The  antenna 
position  and  antenna  current,  i.e.,  the  efficiency  of  the 
transmitting  dipole  were  kept  constant.  The  measured  values 
were  then  plotted  in  a  diagram  over  the  distance  r.  Important 
conclusions  v/ith  respect  to  propagation  and  the  decisive 
factors  can  be  made  from  this  field  strength  profile. 

Another  possibility  which  shall  be  discussed  in  detail 
in  another  section  is  the  measurement  of  an  antenna  diagram. 

Based  on  the  results  of  these  measurements,  the  theore¬ 
tical  team  developed  a  theory  attempting  to  describe  the 
propagation  of  VLF  waves  in  a  dissipative  medium.  In  order 
to  describe  the  method  of  measurement  and  its  results,  the 
theoretical  results  described  in  chapter  4  shall  be 

discussed  in  brief. 

* 

For  a  simpler  mathematical  treatment  of  the  given 
problem,  the  following  assumptions  were  made  for  the  theo¬ 
retical  investigation: 

homogeneous,  isotropic,  unbounded  medium  surrounding 
the  transmitting  antenna  and  the  receiver; 

let  the  transmitting  antenna  be  considered  as  a  magnetic 
dipole; 

let  the  distance  between  the  transmitter  and  receiver 
be  so  large  that  both  may  be  looked  upon  as 
being  point-shaped; 
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lot  us  furthermore  assume  the  decisive  parameters  and 
<fto  be  constant  and  independent  of  the  frequency. 

The  materialization  of  some  requirements,  especially 
that  of  an  unbounded  medium  involves  considerable  difficulties? 
thus,  theoretical  attempts  we re  made,  e.g.,  of  considering 
the  influence  of  the  earth's  surface. 

Since  both,  transmitter  and  receiver,  had  to  be  set 
up  in  adequate  caverns  (shafts,  galleries)  of  the  mine, 
the  influence  of  these  caverns  theoretically  was  estimated 
as  being  spheric'--’  ly  idealized  and  in  most  cases  negligibly 
small  [2}. 

Prom  the  theory  of  a  transmitting  magnetic  dipole, 
the  magnetic  field  strength  was  calculated  which  under 
these  conditions  is  the  only  factor  of  interest.  The 
geometry  of  the  problem  is  placed  in  spherical  coordinates, 
with  the  dipole  axis  coinciding  with  the  z-axis  (cf.  Pig.  2.1). 
With  m  -  pt  *n*I.F  ....  magnetic  moment  of  the  transmitting 

dipole 

n,  I,  F  .  number  of  turns,  antenna  current, 

area  of  the  transmitting  frame, 

rJ~ . angle  between  the  dipole  and  the 

radius  vector  with  respect  to  the 
point  of  measurement, 

r  .  distance  between  transmitter  and 

receiver 

8*  . .  complex  material  constant  containing 

<F,  g,Q,  p,  and  w  =  2m V  (frequency), 

we  thus  have 


•  n* 

o  -18  r 
2m*e  •  cos 


0/1  iS  N 
^  .(~r  +  — ) 

r 


cos 


J  •  h 


H 


=  -m.e 


-18  r 


+ 


(2.1) 

(2.1a) 

I 

(2.2) 

(2.3) 


JET 


—  ~S 
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IhjJ  =  sin  &  .  h^.  (2.3a) 

Prom  =  0  we  obtain  the  symmetry  about  the  z-axis 
(dipole  axis)  in  accordance  with  the  geometry  of  the  trans¬ 
mitting  antenna  which  is  assumed  to  be  circular.  The  field 
lines  thus  all  lie  in  meridional  planes .Hsaoe  the  problem  of 
space  becomes  somewhat  simpler,  the  measurement  of  one 
meridional  plane  is  quite  sufficient  for  describing  m 
detail  the  radiation  picture. 

Only  the  field  strength  value  at  the  respective  point 
of  measurement  which  is  vectorially  composed  of  and  HJJ , 

is  accessible  to  measurement.  The  voltage  U.  ,  which  is 

xnd 

induced  in  the  coil  of  the  receiving  antenna  is  measured 
with  suitably  dimensioned  instruments  (no  overload,  linear 
part  of  the  hysteresis  curve  of  the  ferrite  material); 
it  is  proportional  to  the  field  strength. 


Dind 

-  1  H  1  =  (H 1  +  Hj) 

1/2 

• 

(2.4) 

The  direction  of  the 

field  strength  maximum  a 

ccording  t* 

Pig.  2.2  is  given  by 

l%l 

tan  w  -  . 

15,1 

(2.9) 

Similarly  to  the  Gaussian  positions,  the  measurements  for 

<J*=  0°  and  $■  =  90°  are  of  special  interest. 

From  the  Eq  s. 

(2.1a),  (2.3a)  and 

(2.4a) 

it  follows  that 

fH 

si 

ii 

o 

o 

n 

fa 

(2.6) 

K4 

ii 

vo 

o 

o 

II 

-r 

9s* 

• 

(2.7)) 

Prom  the  relation 

h 

XX 

G  = 

(o 

\ - / 
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conclusions  can  again  be  drawn  as  to  the  decisive  para¬ 
meters,  In  order  to  be  able  to  compare  the  measured  curves 
with  the  theoretical  results,  a  large  array  of  field  strength 
values  was  calculated  for  the  entire  range  of  parameters 
under  consideration,  using  an  electronic  computer  of  the 
Zuse  Z  23  V  type;  the  field  strength  values  were  also 
graphically  represented.  The  curves  shown  in  Pigs. 2. 3  -  2.8 
were  plotted  according  to  a  profile  measurement.  The  constant 
factor  m  (dipole  moment)  was  arbitrarily  taken  as  10^  (in 
accordance  with  the  calculated  and  tabulated  values).  The 
circle  diagrams  in  Pigs.  2.9  -  2.12  correspond  to  the 
well-known  antenna  diagrams  for  certain  distances  and  con¬ 
ductivities.  These  circle  diagrams  were  drawn  for  the  most 
frequently  used  frequency  of  3  kcps. 

2.1.1.  Methods  of  measurement 

Profile  measurement 

This  very  simple  possibility  of  measuring  a  VLP  field 
has  already  been  described  at  the  beginning.  Prom  it  re¬ 
sults  the  magnitude  of  the  decrease  in  field  strength  with 
increasing  distance  r. 

Antenna  diagram 

An  antenna  diagram  is  measured  in  the  well-known 
manner  by  finding  the  field  strength  values  at  a  constant 
distance  around  the  transmitting  antenna*,  if  necessary 
also  their  direction  is  measured.  As  such  measurements 
are  impossible  in  a  mine,  the  transmitting  antenna  was 
turned  instead,  and  the  receiving  antenna  was  left  at  its 
place.  The  rotation  of  the  transmitting  antenna  may  be 
seen  from  the  angle  $  in  Pig.  2.2.  The  direction-finder 
antenna  PA  5  (for  details  see  Chapter  l)  was  developed  for 
the  purpose  of  measuring  the  field  strength  values  at  the 
point  of  measurement,  as  well  as  for  a  precise  angular  measure- 
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raent  of  the  field  direction.  The  ferrite  rod  is  mounted  on 
a  turn-table  so  as  to  rotate  perpendicularly  to  its  axis. 
The  turn-table  in  turn  swivels  in  any  direction;  for  field 
strength  measurements  it  is  adjusted  so  that  its  plane 
contains  the  field  strength  maximum.  The  latter  determines 
the  angle  ^  (see  Fig.  2.2).  In  order  to  reduce  the  errors 
in  reading  and  adjusting  the  angles,  another  circle  diagram 
is  recorded  with  the  use  of  tfc  receiving  antenna:  The 
ferrite  rod  is  rotated  through  10°  each  time'etarting  from 
the  maximum  direction  ( =  0°),  and  the  respective  voltage 
values  are  read;  between  the  latter  and  the  field  strength 
maximum  there  exists  the  following  relation: 


ur  =  uind*  cosr 


(2.9) 


U.  ,  is  determined  for  every  $  by  a  balancing  calculation 
ind 

(method  of  least  squares)  and  from  it  the  antenna  characte¬ 
ristics  are  concluded  by  another  balancing  calculation. 

This  antenna  characteristic  having  the  form  of  a 
circle  diagram  (Figs.  2.9  -  2.12)  also  contains  the  radial 
field  strength  h^  for  ^  =  0°  and  for  =  90°  the  tangential 
field  strength  h^  .  The  ratio  G  =  h^:  h&  is  a  near  field  - 
-  far  field  measure  depending  considerably  on  the  conducti¬ 
vity  [3],  In  the  extreme  near  field  (r  very  small). 

G  =  2  (cf.  Eqs.  (2.1)  and  (2.3)).  The  relation  shifts  toward 
smaller  G  values  as  the  distance  increases.  From  G  <  1  on¬ 


ward  we  may  speak  of  a  far  field.  In  connection  with  it,  the 
maximum  field  strength  of  reception  in  its  vicinity  can  be 
expected  at  an  angle  of  <^=  0°,  at  greater  distances  (the 
boundary  varies  depending  on  the  conductivity  and  frequency) 
the  principal  radiation  direction  is  ^  90°,  a  minimum 
occurring  at  J^=  0°.  For  G  =  1,  the  radiation  characteristic 
is  spherical. 


Measurement  of  angle 


In  the  method  described  so  far,  the  propagation  para- 
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meters  (the  most  important  one  being  the  conductivity  6  )  in 
general  are  determined  only  as  a  value  integrated  over  the 
entire  range  of  measurement.  For  various  purposes  of  applic¬ 
ation  it  is  desirable  that  disturbances  even  smaller  than  that 
be  detected  which  integration  over  the  entire  region  would 
be  too  small  as  to  affect  the  result  noticeably.  This  is 
valid  especially  in  those  cases  where  the  attenuating  effect 
of  the  dissipative  medium  at  the  distance  in  question  is  not 
noticed  (cf.  also  Figs.  2.3  -  2.8,  linear  section  of  the 
curves)  because  of  the  low  conductivity. 

Disturbances  o.s  small  as  that  do  not  affect  the  amplitude 
to  a  measurable  extent,  a  deviation  of  the  angular  position 
of  the  maximum,  however,  may  be  expected  j_ 4 J  .  From  Eq-.  .(2.5) 
(cf.  Fig.  2.2)  we  obtain: 


h  ft  *  sm  'S-  0 

,  v-  tan 

tan  u/  =  - - g-  =  — r — 

T  h  •  cos  a>-  G 


(.2.10) 


The  dependence  of  the  angle  y  on  for  different  values  of 
G  is  shown  in  Fig.  2,13. 

The  direction  of  the  maximum  is  determined  not  at  a 
weak  maximum,  but  at  the  clearly  adjustable  minimum  per¬ 
pendicular  to  it.  The  direction-finder  antenna  FA  5  con¬ 
structed  for  this  purpose,  permits  an  accuracy  of  measure¬ 
ment  <  1°, 


2.1.2.  Modeling  tests  [5j 

Within  the  framework  of  experimental  studies,  modeling  tests 
were  also  taken  into  consideration  [3,  4,  5]  .  it  would  be 
highly  desirous  that  inhomogeneities  in  the  rock  should  be 
avoided  and  that  the  conditions  of  laboratory  tests  should  be 
accurately  expressed  by  theory.  This  is  of  great  interest 
especially  in  those  cases  where  the  earth  -  air  interface  is 
being  considered,  which  in  most  mines  is  not  sufficiently  plane, 
especially  not  in  rhe  Tyrolean  mines. 
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Another  advantage  would  be  the  independence  of  accessible 
galleries  that  are  suitable  for  our  measurement;  the  choice 
being  highly  restricted  above  all  by  the  latter  condition. 

For  a  modeling  test  it  is  by  no  means  sufficient  to 
change  all  longitudinal  dimensions  (wavelength  and  with  it  the 
frequency),  but  the  conductivity  and  the  dielectric  constant 
must  also  be  varied  [23j  « 

Based  on  the  wave  equation  we  divide  all  quantities  into 
factors  with  and  without  dimensions: 

•* 

A’E 1  -  (r2p£m2  -  itor2Sp)E'  =  0  (2.1l) 

(dashed  quantities  being  dimensionless). 

This  dimensionless  equation  is  not  changed  any  further 
by  a  gauge  transformation,  as  long  as  the  individual  terms  are 
kept  constant,  which  themselves  are  not  dimensionless.  Thus 
we  obtain  the  "modeling  condition"  (subscript  m  ...  model  size): 


real  part:  r2p£w2  =  r2  p  £  u)2 
^  r  mm  m  m 


imaginary  part:r2p6u>  =  r2  p  o' 

mm  m  m 


(2.12) 


Since  the  electrical  quantities  of  the  media  in  question 
differ  widely,  we  may  distinguish  two  cases: 

Case  1: 


x 

Sw 


»  1 


Case  2: 


«  1 


(imaginary  part  of  the  wave  equation 
prevails,  high  conductivity,  e.g., 
earth,  ionosphere).  (2.13) 

(Real  part  prevails,  vanishing  conduc¬ 
tivity,  e.ge;air)«  (2=14) 


In  both  cases  only  one  oondition  (2.12)  must  be  satisfied, 
which  is  an  essential  simplification. 

For  an  earth  -  air  -  ionosphere  model,  the  meaningful 
model  factors  have  an  order  of  10b .  From  it  results  a  model 
conductivity  which  can  be  represented  by  metal  only,  "Underground" 
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secondary  effects  (skin  effect)  occur  at  the  model  frequencies 
(several  Mcps),  rendering  such  a  measurement  impossible. 

Reducing  the  model  factor  such  that  the  earth  can  be  re¬ 
presented  by  acidified  water  makes  the  dimensions  un¬ 
suitable  . 

Another  field  of  application  would  be  measuring  the 
antenna  directivity  patterns  under  clearly  defined  con¬ 
ditions.  In  this  case,  the  model  factors  are  approximately 
2  3 

10  to  10^,  thus  involving  conditions  that  can  be  materialized. 

Furthermore  the  possibility  was  studied  of  making 
modeling  tests  on  solid  rock  in  the  mine  in  order  to  be 
able  to  conduct  measurements  over  great  distances  also  at 
restricted  regional  dimensions.  At  constant  conductivity, 
only  length  and  frequency  must  be  changed. 

Case  1:  The  simultaneous  fulfilment  of  Eqs.  (2.12)  and 
(2.13)  (imaginary  part)  restricts  the  number  of 
possibilities.  The  condition 

•—  =  —  »  1  (2.13) 

puts  an  upper  limit  of  frequency  of  approximately  10  kc/sec 
at  &  =  10  ^  mhos/n.  This,  however,  lies  in  the  order  of 
magnitude  of  the  actual  frequency  of  measurement.  A  model 
of  this  type  is  possible  only  at  higher  conductivities. 

Case  2:  The  condition  (2.12)  (real  part)  is  independent 
T  of  the  conductivity. 

According  to  the  requirement  of  (2.I4),  the  latter  con¬ 
dition  is  restricted  to  conductivity  values  below  10~^  mhos/m 
and  frequencies  above  1  Mc/sec. 

For  the  region  that  lies  between  case  1  and  case  2, 
which  is  quite  frequent  in  practical  application,  the  fulfillment 
of  both  conditions  (2.12)  is  necessary?  this  is  possible 
only  when  length,  frequency  and  conductivity  are  modeled 
s imul taneous ly . 
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The  only  possibility  of  conducting  laboratory  tests  at 
reasonable  modeling  dimensions  (having  the  order  of  meters) 
is  the  measurement  of  antenna  directivity  patterns,  slotting 
such  diagrams  proved  to  be  comparatively  easy  in  practice, 
especially  since  in  numerous  cases  a  distance  of  less  than 
1  km  was  measured  so  that  the  expenses  for  setting  up  bhe 
instruments  and  for  the  modeling  tests  themselves  did  not 
appear  to  be  justified. 

2.2.  Measurements  in  the  minp 

The  measurements  below  ground  were  made  in  several  mines 
offering  various  types  of  rock.  Thus  it  w as  possible  to 
obtain  series  of  measurements  for  a  large  range  of  decisive 
parameters  which  were  then  compared  with  the  theoretical 
results  obtained  for  wave  propagation  in  solid  media  ij  , [.2J  . 

The  conductivity  of  the  surrounding  rock  varies  betv/een 
approximately  10  ^  mhos/m  in  the  mine  of  Konrad  I  (i-rOn  ore) 
to  10-8  mhos/m  in  the  Hansa  III  mine  (dry  salt). 

The  sites  of  measurement  v/ere  chosen  so  that  the  conditions 
assumed  for  theoretical  calculation  were  fulfilled  as  far 
as  possible.  Large  iron  units  such  as  rails,  pipes, 
mining  railways  or  iron  formwork  had  to  be  avoided.  In  order 
to  meet  the  requirement  of  a  homogeneous  medium,  there  should 
be  no  dislocations,  cavities  or  ore  inclusions  in  the  region 
of  measurement. 

The  distances  between  the  individual  points  of  measure¬ 
ment  were  read  from  mine  charts  and  the  angles  were  measured  ' 
with  a  geological  compass.  The  resulting  errors  of  measure¬ 
ment  (inaccurate  charts,  magnetic  disturbance  of  the  compass) 
are  partly  corrected  by  balancing  calculations. 

The  necessary  communication  between  the  transmitter  and 
the  sites  of  reception  is  made  possible  by  a  field  telephone. 
The  effect  of  the  connecting  wire  on  the  measurement  can 
be  neglected  according  to  various  experiments.  In  cases 
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where  a  telephone  communication  was  impossible  because  of 
the  great  distances,  the  measurements  were  made  in  accordance 
with  an  accurate  timetable.  The  distance  over  which  measure¬ 
ments  were  conducted  varied  considerably,  depending  on  the 
conductivity  and  the  possibilities  of  the  respective  mine. 

In  many  cases,  measurements  could  be  made  over  distances  of 
800  m,  occasionally,  however,  only  20  to  30  m.  In  the  "home 
mine"  of  St.  Gertraudi,  Tyrol,  measurements  over  a  distance 
of  7000  m  with  a  stationary  transmitter  erected  for  this 
purpose  were  possible,  part  of  the  distance,  however, 
being  above  ground. 

2.2.1.  Mine  of  Lafatsch,  Tyrol 
(Mining  charts  see  [2]) 

This  mine  which  at  present  has  been  abandoned  contains 
a  large  deposit  of  lead  and  zinc  ores  in  the  form  of  a  steep 
layer,  entered  at  several  levels.  The  transmitter  was  first 
set  up  at  the  lowest,  the  sixth  level  having  the  largest 
system  of  galleries,  where  the  majority  of  measurements  were 
conducted.  The  ore  is  mainly  located  above  the  region  of 
measurement . 

The  conductivity  determined  on  the  basis  of  these  measure¬ 
ments  by  a  balancing  calculation  is  approximately  5«10  ^mhos/m. 

For  a  second  series  of  measurements,  the  transmitter 
was  set  up  at  the  top  end  of  the  ore  layer  on  level  2,  the 
field  strength  measurements  being  conducted  along  the  ore 
deposit  down  to  level  6.  As  could  be  expected  from  the 
theoretical  results,  the  integrated  conductivity  was  some¬ 
what  higher,  namely  10  ^  mhos/m.  Both  curves  together  with 
the  theoretical  representation  belonging  to  this  conductivity 
value  are  shown  in  Fig.  2.14.  It  is  quite  evident  that  the 
ore  deposit  swallows  more  energy  because  of  its  higher 
conductivity  (curve  l)  so  that  the  curve  tr^nd  on  the  whole 
is  somewhat  steeper  for  measurements  along  the  ore  deposit. 
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Although  the  conductivity  during  the  measurements  at 
Lafatsch  was  found  to  differ,  the  maximum  distance  of  measure¬ 
ment  at  the  given  frequency  and  material  constants  was  too 
small  as  to  yield  a  satisafactory  accuracy  of  conductivity 
values . 


2.2.2.  Mine  at  Schwaz ,  Tyrol 

(Mining  charts  see  [2]) 

In  contrast  to  the  Lafatsch  measurements,  the  measurements 
here  oovered  distances  of  up  to  approximately  1300  m.  For 
this  distance,  the  previous  transmitter  was  no  longer  suited. 
Another  transmitting  antenna  of  about  16  ma  and  35  turns 
(for  details  of  SA  VII  see  [2]  ),  was  set  up  for  covering 
this  distance.  The  measured  values  are  plotted  in  Fig.  2.15* 
Curve  3  was  measured  with  the  small  transmitting  antenna 
(SA  III),  curve  2,  however,  shows  the  values  obtained  with 
the  large  frame.  The  measurements  of  curve  1  were  made  through 
different  types  of  rock  (slate  wedge,  else  dolomite).  The 
mean  conductivity  calculated  from  it  is  10-^  mhos/m. 

The  measurements  made  so  far  were  almost  pure  near  field 
measurements  owing  to  the  small  conductivity  and  distances. 

They  were  consistent  with  the  theory  within  the  measuring 
accuracy,  although  the  compared  values  were  not  absolute, 
a  balancing  calculation  was  made  in  order  to  make  the 
measured  attenuation  consistent  with  the  theory. 

2.2.3*  Konrad  I  mine  at  Salzgitter 

The  geological  structure  of  this  region  is  characterized 
by  a  varied,  flat  stratification  formed  by  repeated  sea  floods. 
The  concentration  of  this  oolitic  ore  is  approximately  30$, 
the  average  size  of  the  individual  brown  iron  grains  is  0.5  mm. 
The  depth  of  the  mine  is  approximately  1000  m. 

Sites  suitable  for  the  measurements  were  difficult  to 
find  in  the  mine,  they  had  to  be  free  of  iron  formwork,  rails, 
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pipes  and  power  cables.  An  18  m  and  a  25  a  distance  at  the 
1000  m  level  appeared  to  be  suited  best.  Although  the  distances 
were  short,  characteristic  antenna  diagrams  could  still  be 
measured,  as  the  conductivity  was  very  high,  namely  2 .10~^mhos/m. 

The  measurements  in  Konrad  I  showed  that  the  voltage  read 
at  an  angle  of  ^  =  90°  (cf.  chapter  2.1.1.)  does  not  always 
become  zero  as  had  been  expected  according  to  the  theory  of 
the  cosine  law.  An  explanation  although  just  a  qualitative 
one,  would  be  the  high  inhomogeneity  of  the  field,  an 
assumption  which  seems  quite  probable  [3J  at  the  small  distance 
and  the  high  conductivity.  A  detailed  examination  of  this  effect 
will  be  made  at  higher  frequencies. 

An  evaluation  of  the  antenna  diagrams  showed  that  the 
transition  from  near  field  to  far  field  at  2*10  ^  mhos/m 
occurs  already  at  25  m  (G  ~  l),  whereas  at  Gertraudi 
(approximately  10*^  mhos/m)  the  measurements  were  still  made 
in  the  near  field  even  at  800  m  (G  >  l).  In  this  case  of 
extremely  high  conductivity  the  strong  shift  in  near  field  - 
far  field  transition  becomes  quite  evident.  Under  these 
circumstances,  effects  occur  at  minimum  distances  which 
elsewhere  set  in  at  greater  distances,  being  immeasurable 
because  of  the  small  power  of  the  transmitter. 

2.2.4.  Potassium  mine  Hansa  III,  Empelde  near  Hannover 

Uhereas  in  Konrad  I  the  electrical  conductivity  of  rock 

on  average  was  by  three  powers  higher  than  in  the  mines  studied 

so  far,  the  conductivity  of  the  medium  used  in  Hansa  III  was 

—8 

so  small  (approximately  10  mhos/m)  that  it  may  be  looked 
upon  as  being  almost  an  insulator.  The  measurements  were 
conducted  in  a  large,  dry  salt  stock  at  the  1000  m  level. 

Apparently  favorable  conditions  of  propagation  ana 
measurement,  however,  could  not  be  used  much  because  of  the 
disturbing  influences.  Owing  to  the  weak  attenuation  in  the 
salt  body,  disturbances  caused  by  the  mining  railway  at  higher 
levels  were  noted  stronalv.  High  static  charges  caused  bv 
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pipes  and  power  cables.  An  16  m  and  a  25  m  distance  at  the 
1000  m  level  appeared  to  be  suited  best.  Although  the  distances 
were  short,  characteristic  antenna  diagrams  could  still  be 
measured,  as  the  conductivity  was  very  high,  namely  2.10~'*'mhos/m. 

The  measurements  in  Konrad  I  showed  that  the  voltage  read 
at  an  angle  of  ^  -  90°  (cf.  chapter  2.1.1.)  does  not  always 
become  zero  as  had  been  expected  according  to  the  theory  of 
the  cosine  law.  An  explanation  although  just  a  qualitative 
one,  would  be  the  high  inhomogeneity  of  the  field,  an 
assumption  which  seems  quite  probable  [3J  at  the  small  distance 
and  the  high  conductivity.  A  detailed  examination  of  this  effect 
will  be  made  at  higher  frequencies. 

An  evaluation  of  the  antenna  diagrams  showed  that  the 
transition  from  near  field  to  far  field  at  2»10~‘1'  mhos/m 
occurs  already  at  25  m  (G  ~  l),  whereas  at  Gertraudi 
(approximately  10"^  mhos/m)  th  >  measurements  were  still  made 
in  the  near  field  even  at  800  m  (G  >  l).  In  this  case  of 
extremely  high  conductivity  the  strong  shift  in  near  field  - 
far  field  transition  becomes  quite  evident.  Under  these 
circumstances,  effects  occur  at  minimum  distances  which 
elsewhere  set  in  at  greater  di stances  being  immeasurable 
because  of  the  small  power  of  the  transmitter. 
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V/hereas  in  Konrad  I  the  electrical  conductivity  of  rock 
on  average  was  by  three  powers  higher  than  in  the  mines  studied 
so  far,  the  conductivity  of  the  medium  used  in  Hansa  III  was 
so  small  (approximately  10  mhos/m)  that  it  may  be  looked 
upon  as  being  almost  an  insulator.  The  measurements  were 
conducted  in  a  large,  ary  salt  stock  at  the  1000  m  level. 

Apparently  favorable  conditions  of  propagation  and 
measurement,  however,  could  not  be  used  much  because  of  the 
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salt  body,  disturbances  caused  by  the  mining  railway  at  higher 
levels  were  noted  strongly.  High  static  charges  caused  by 
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salt  yielded  considerable  irregular  fluctuations  in  indication. 
The  GBR  transmitter  which  above  ground  was  very  strong,  could 
not  be  received,  since  the  salt-containing  ground  water  layer 
at  a  depth  of  180  m  screens  off  all  signals  from  above. 

Yet,  an  antenna  diagram  was  measured  whose  ratio 
G  -  h  :  h^  =  2  leads  to  the  assumption  of  a  very  low  con¬ 
ductivity. 

2,2.5*  Ore  mines  of  Fusseberg  and  Georg 
(Sieger land,  Germany) 

From  an  electric  viewpoint,  the  conditions  in  these 
two  mines  were  similar  to  those  at  Sbhwaz  or  Gertraudi.  The 
iron  ore  (siderite)  mined  there  occurs  in  the  form  of  a 
steep  edge  heave  approximately  10  m  thick.  The  deposit  of 
Fusseberg  is  straight,  whereas  in  the  Georg  mine  it  has 
the  shape  of  a  horseshoe. 

The  sites  of  measurement  were  chosen  so  that  the  line 
connecting  transmitter  and  receiver  ran  mainly  through  ore. 

The  influence  of  the  inhomogeneous  medium  on  the  propagation 
of  VLF  waves  could  thus  be  studied.  The  measured  values 
showed  on  the  whole  the  same  trend  as  those  measured  at 
Schwaz,  thus  suggesting  equal  conductivity.  Direct  current 
measurements  which  were  made  parallel  [jJ  yielded  the  same 
conductivity  for  ore  and  dead  rock,  thus  explaining  the  small 
differences  on  measuring  in  front  of  and  behind  the  ore  layer. 

The  antenna  diagrams  shov/ed  pure  minima  at  ^/  =  90°  and 
were  well  consistent  with  the  theoretical  values  obtained  for 
this  conductivity  value. 

2.2.6.  Goal  mines  Oran.je  Nassau,  Heerlen  (Netherlands) 

The  Nassau  III  mine  offered  the  unique  possibility  of 
making  measurements  in  coal.  They  were  made  in  the  direction 
of  the  ooal  seam  which  is  about  2  m  thick.  For  this  purpose, 
the  mine  management  provided  a  surrounding  for  transmitter  and 
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receiver  that  was  free  of  iron,  short  drifts  were  lined  with 
wood  formwork.  Measurements  were  made  at  the  angles  ft  =  0* 
and  /3~  =  90°  and  at  a  distance  of  300  m  and  500  m,  respectively. 
Unexpectedly,  the  conductivity  measurements  with  direct  current 
showed  coal  to  have  a  poor  conductivity  and  the  rock  below 
was  found  to  have  a  comparatively  high  conductivity.  Hence, 
propagation  proceeded  in  the  direction  of  a  thin  and  weakly 
conductive  layer  which  was  bounded  by  two  media  of  good  con¬ 
ductivity.  Such  a  "sandwich"  structure  which  in  principle  is 
typical  of  any  coal  pit,  on  account  of  the  results  seems  to 
favor  the  propagation  of  VLF  waves.  Various  factors  apparently 
affect  the  propagation:  The  relation  between  the  seam  thick¬ 
ness  and  the  measuring  distance  in  connection  with  the  wave¬ 
length  which  again  depends  on  the  conductivity  of  the  surrounding 
medium,  is  an  essential  factor.  Also  the  antenna  alignment 
(angle  /ib)  has  a  different  effect  as  compared  to  air  or  homo¬ 
geneous  rock. 

At  a  distance  of  500  m,  G  was  found  to  have  the  value 
0.22  which  corresponds  to  a  conductivity  of  about  l»10~^mhos/m 
at  a  frequency  of  3  kepa.  A  comparison  with  the  d-c  conductivity 
however,  cannot  be  made,  3ince  the  propagation  in  the  strongly 
heterogeneous  medium  cannot  be  compared  with  the  conditions 
of  our  theory.  It  was  not  possible  to  study  the  variety  of 
effects  on  the  spot  7*  1  tii  *  a*  e  required  accuracy.  The  prelimi¬ 
nary  results  show,  however,  that  a  detailed  examination  of 
this  alignment  would  yield  highly  interesting  results, 

2.2.7.  Mine  at  Blelberg,  Karnten  (Austria)  [15J 

In  this  mine,  measurements  were  mainly  made  in  Yfetter- 
stein  limestone  containing  ore  inclusions  (20  to  30  m)  in 
the  form  of  lead  glance  and  zinc  blende  which  compared  to  the 
range  of  measurement  (400  to  800  m)  were  rather  smell.  The 
inhomogeneities  therefore  had  only  a  slight  effect  on  the 
results  of  measurement.  Fig.  2.16  shows  such  a  profile  measure- 
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ment  (curve  l).  The  dashed  curve  gives  the  theoretical  decrease 
in  field  strength  resulting  from  the  mean  value  of  conductivi¬ 
ty,  approximately  10*"^  mhos/m.  At  about  400  m,  the  transition 
to  rock  of  somewhat  higher  conductivity  is  clearly  recognizable. 

Curves  2  and  3  give  the  measurements  of  the  GBR  trans- 
mitter  (l6  kcps).  The  horizontal  field  strength  decrease  was 
plotted  for  the  distance  from  the  shaft  mouth  (curve  2)  and 
ir.  2a  the  cover  is  also  plotted  (shortest  distance  from  the 
earth's  surface).  Curve  3  gives  the  decrease  in  vertical 
direction.  The  distance  here  is  nob  the  transmitter  -  receiver 
range  which  ha3  the  order  of  about  1000  km,  but  the  respective 
distance  from  the  earth's  surface.  Hence,  the  much  weaker 
attenuation  of  the  GBR  transmitter  as  compared  to  our  own 
transmitter  can  also  be  explained,  since  in  one  case  the 
measurement  is  made  in  the  far  field  (approximately  plane 
waves),  whereas  in  the  other  case  the  geometry  of  the  field 
contributes  considerably  to  the  attenuation  (decrease  of 
the  field  with  l/r^  in  addition  to  absorption). 

2 .2 .8  .  Mine  at  St.  Gertraudi.  Tyrol 

(mining  charts  see  [2]  and  [3J  ) 

At  Gertraudi,  a  straight  gallery  300  m  long  was 
available  as  a  connection  between  the  GroBkogei  mine  and  the 
Kleinkogel  mine,  which  is  highly  suited  as  it  passes  through 
largely  homogeneous  dead  rock;  its  elongation  of  another  500  m 
in  the  abandoned  Kleinkogel  mine  is  again  highly  suitable 
for  the  measurements.  This  is  the  reason  why  preliminary 
experiments  and  measurements  were  made  in  this  mine  from 
the  very  beginning  of  our  project  work,  besides  the  measure¬ 
ments  made  at  Schnmz, 

Let  me  mention  some  of  the  numerous  measurements  and 
experiments:  Preliminary  experiments  with  adapted  transistorized 
portables;  measurements  of  other  long-wr.ve  transmitters; 
primary  experiments  v/ith  our  own  transmitter  at  higher  fre- 
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quencies  (about  100  kcps);  profile  measurements  with  specially- 
designed.  transmitting  antennas;  experiments  and  measurements 
V7ith  various  receiving  antennas. 

On  the  preparatory  measurements  discussed  herein,  a  theory 
of  wave  propagation  was  based  which  is  adapted  to  the 
simplified  conditions  of  measurement  [2l  as  has  been  mentioned 
in  2.1.  Prom  the  profile  measurements,  the  conductivity  was 
calculated  and  found  to  be  of  the  order  of  10  ^  mhos/m,  in 
some  series  of  measurements,  however,  only  the  conductivities 
above  10~^  mhos/m  at  smaller  distances  could  be  excluded.  A 
comparison  with  the  theory  at  first  could  only  be  based  on 
the  attenuation  relation  owing  to  the  absence  of  accurate 
calibration  factors  for  transmitting  antennas  and  receiving 
antennas,  i.e.,  the  slopes  of  the  measured  curves  were  compared 
with  the  respective  theoretical  curves  by  means  of  a  balancing 
calculation  (method  of  least  squares);  the  agreement  with 
theory  was  good  so  that  the  adapted  theory  of  the  transmitting 
dipole  for  the  near  field  may  be  assumed  to  be  correct  [ 2 J  . 

The  antenna  diagrams  measured  in  the  direct  neighborhood 
of  the  transmitting  antenna  also  showed  good  agreement  with 
the  theory  for  transmitting  dipoles. 

Prom  the  electrical  and  geological  viewpoints,  the  con¬ 
ditions  of  measurement  in  the  mines  of  Schwas  and  St.  Gertraudi 
are  approximately  equal,  hence  also  the  similar  results  of 
measurements  (cf.  chapter  2.2.2.). 

Prom  the  requirement  for  covering  larger  distances  re¬ 
sulted  the  necessity  of  building  always  larger  transmitters 
and  of  making  the  receiving  antennas  more  selective  and 
sensitive.  By  means  of  a  stationary  1  kw  amplifier  set  up 
in  the  mine  laboratory  [33  (details  see  chapter  l),  distances 
of  up  to  7  kin  could  be  covered  by  means  of  the  transmitting 
antenna  SA  VIII  (dimensions  40  x  40  m.)eet  up  especially  for 
this  purpose.  As  a  distance  Of  the  above  length  could  not  be 
found  within  the  mine,  measurements  at  the  remotest  points  had 
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to  be  made  above  ground.  Owing  to  the  much  higher  noise 
level  (no  screening  effect  by  rock  layers),  the  sensitivity 
of  the  receiving  antenna  could  not  be  used  fully  so  that  the 
corresponding  measurements  below  ground  would  permit  even 
larger  distances  [3I  . 

An  improved  version  of  the  SA  VII  transmitting  antenna 
was  used  for  plotting  antenna  diagrams  at  frequencies  of  3 
and  11  kcps.  As  the  distances  covered  with  the  SA  VII  antenna 
are  shorter  than  those  of  SA  VIII,  diagrams  were  made  up 
to  a  distance  of  no  more  than  820  m.  The  diagram  of  Fig.  2.l6a 
clearly  shows  that  the  measurement  was  made  in  the  near 
field  also  at  greater  distances,  thus  being  in  good  agree¬ 
ment  with  the  theory  for  the  existing  conductivity  values  |_4_j. 
In  this  connection,  the  change  in  angular  position  ^ of  the 
receiving  antenna  was  also  studied.  The  corresponding  points 
are  plotted  in  Fig.  2.13.  Here,  agreement  with  the  theory 
is  not  quite  as  good,  as  small  inhomogeneities  affect  the 
angular  position,  whereas  the  influences  on  the  amplitude 
cannot  be  measured  (cf.  chapter  2.l). 

2.3.  Comparison  between  experiment  and  theory 

The  great  number  of  measurements  which  were  made  under 
different  geological  conditions  were  almost  all  confined 
to  the  near  field,  with  the  corresponding  distance  of 
measurement  changing  considerably  in  dependence  on  the 
conductivity.  It  has  already  been  said  in  chapter  2.1  that 
propagation  at  very  small  distances  is  quite  independent 
of  external  influences.  The  results  of  measurement  therefore 
are  in  good  agreement  with  the  theory,  since  the  assumptions 
on  which  the  theory  is  based,  are  actually  verified  (measure¬ 
ments  of  [l]  and  partly  also  of  j  2 J  ), 

For  greater  distances  (measurement  still  in  the  near 
field),  the  influence  of  the  conductivity  £  is  already 
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noticeable;  it  is  therefore  used  for  determining  the  value 
of  &  .  In  this  case,  the  results  of  measurement  are  satis¬ 
factorily  explained  by  the  theory  (measurements  of  chapters 
2.2.1,  2.2.2,  2.2.4,  2.2.5). 

In  few  cases  the  conditions  of  measurement  did  not  at 
all  correspond  to  the  theoretical  conditions  (measurement  in 
the  coal  seam,  chapter  2.2,6),  thus  making  a  comparison 
impossible . 

The  measurement  in  the  mine  Konrad  I  (chapter  2.2.5)  was 
made  in  a  medium  of  extremely  high  conductivity  (  &  =  2  •10*'^' 

mhos/m),  yielding  for  the  first  time  a  result  that  can  be 
theoretically  interpreted  in  terms  of  a  far  field  measure¬ 
ment  (G  ~  0.2) . 

The  most  recent  measurement  over  large  distances  at 
St.  Gertraudi  (chapter  2.2.8)  offered  the  possibility  of 
studying  also  the  effects  of  reflection  from  the  earth's 
surface  and  surface  waves  (cf.  Pig.  2.17).  A  comparison 
with  the  theory  of  [lol  (reflection  from  the  earth  -  atmosphere 
interface)  is  only  conditionally  possible,  as  the  earth's 
surface  in  the  region  of  the  St.  Gertraudi  mine  is  by  no 
means  plane,  as  assumed  for  the  present  calculation.  Yet, 
propagation  in  accordance  with  the  refined  theory  is 
quite  probable.  Slimming  up  the  results  we  may  say  that  the 
propagation  of  VLF  waves  in  the  near  field  (G  >  l)  on  the 
whole  is  described  by  the  theory  developed  in  [lj  and  jjfj  . 

For  larger  distances,  this  theory  probably  is  too  simple, 
a  correction  according  to  [loj  (reflection  from  the  earth's 
surface)  yields  better  agreement.  An  influence  of  the  medium 
in  the  form  of  refz  action  or  diffraction  of  VLF  waves  at  the 
existing  power  and  frequency  of  transmitter  could  not  be 
observed  to  a  useful  extent.  At  present,  comprehensive 
experiments  are  made  at  approximately  100  kcps  fox*  studying 
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2.4*  Bearing  of  a  dipole  probe 

(an  example  of  practical  application  [ 20~]  ) 

From  the  determined  experimental  and  theoretical 
results,  the  possibility  of  taking  a  bearing  of  a  transmitting 
magnetic  dipole  in  a  solid  medium  with  high  accuracy  was 
derived.  From  the  expressions  (2.1),  (2.2),  and  (2.3) 

(chapter  2.1)  it  follows  that  all  field  lines  lie  in  me¬ 
ridional  planes,  i.e.,  every  tangent  to  a  field  line 
crosses  the  dipole  axis  (cf.  Fig.  2.18).  This  tangential 
direction  is  exactly  the  direction  at  which  maximum  voltage 
is  induced  in  the  ferrite-rod  direction-finding  antenna. 

In  general,  this  direction  does  not  directly  point  to  the 
transmitting  dipole,  as  can  be  seen  in  the  field  line 
diagram  of  Fig.  2.20a.  If  the  magnetic  dipole  is  displaced 
along  its  axis  (position  B),  the  direction-finder  antenna 
shows  a  new  maximum  direction  b  at  the  same  point  of 
reception.  The  two  space  directions  a,b  define  one  me¬ 
ridional  plane  (two  intersecting  lines). 

For  determining  another  meridional  plane,  two  more 
field  directions  are  measured  from  the  point  II  (cf.  Fig.  2.19) 
in  the  same  manner  as  above  (c,d  correspond  to  the  positionc 
C,D).  The  positions  C,  D  may  coincide  with  A,  B,  but  can 
be  chosen  arbitrarily  in  accordance  with  the  respective 
conditions.  The  line  of  intersection  of  these  two  meridio¬ 
nal  planes  is  identical  v/ith  the  dipole  axis. 

The  position  of  the  basis  line  I-II  with  respect  to  the 
dipole  axis  must  be  wind-tipped,  a  parallel  position  or 
intersection  of  the  cwo  lines  would  render  the  problem  in¬ 
definitive  . 

The  accurate  location  of  $ie  dipole  itself  after  having 
found  its  axial  direction  by  means  of  the  field  line  diagram 
(Fig.  2,20  a)  is  now  possible.  The  direction  of  field  lines 
is  given  by  the  expression: 
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with  the  expressions  for  and  (expressions  (2.l)  and 
(2.3),  chapter  2.l)  being  approximated  by  the  first  term, 
in  accordance  with  the  problem  (measurement  over  a  compara¬ 
tively  short  distance). 

Integration  of  the  equation  (2.1?)  yields 

r  =  a  •  sin^  $■  (2.16) 

for  the  field  lines,,  v/ith  a  being  the  parameter  of  field 
lines . 

Similar  to  the  antenna  diagram  measurements  of  chapter 
2.1.1,  the  displacement  of  the  transmitting  dipole  (Pig. 

2.20  d)  with  the  corresponding  directions  a,h  here  is  re¬ 
presented  as  a  shift  of  the  point  A,B,  thus  permitting 
a  full-scale  representation  of  the  direction  of  reception  in 
one  diagram  of  field  lines  (Pig.  2.20  a).  As  the  angular 
position  of  the  field  lines  with  respect  to  one  radius 
vector  is  always  equal  and  independent  of  the  radius,  the 
scale  may  be  chosen  arbitrarily.  For  a  unique  determination 
of  position,  always  two  angular  positions  must  be  taken  into 
consideration. 

The  method  described  above,  is  suited,  e.g.,  for  lo¬ 
cating  miscarried  boreholes  in  mines.  Because  of  the  narrow 
hole,  the  magnetic  dipole  for  this  purpose  usually  has  the 
shape  of  a  long  and  thin  coil  with  an  iron  core  placed  at 
various  depths  (according  to  the  positions  A...D). 

In  a  measurement  based  on  this  method,  the  position  of 
a  borehole  (54  nun  diameter)  which  was  miscarred  by  1.6  m  in 
the  copper  mine  of  I/Iitterborg,  Salzburg,  could  be  determined 
with  an  accuracy  of  10  cm  f  20l . 
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3.  Studies  on  the  electrical  conductivity  of  rock 

The  studies  on  the  propagation  of  VLF  waves  also  included 
measurements  of  the  conductivity  and  dielectric  constant  of 
rock.  The  conductivity  measurements  on  solid  rock  were  made 
with  direct  current  and  alternating  current,  whereas  those 
on  rock  samples  were  made  with  alternating  current  only. 

3.1.  Direct  current  measurements 


The  conductivity  measurements  with  direct  current 
followed  the  principles  of  the  four-electrode  configuration. 
The  theoretical  fundamentals  and  the  measuring  equipment 
used  for  this  purpose  are  described  in  [  $~j  ,  p.  1-23. 

The  principal  aims  of  direct-current  conductivity  measure¬ 
ments  were  the  following: 

(a)  to  study  how  far  rock  is  uniform  from  an  electrical 
point  of  view, 

(b)  to  derive  the  conditions  under  which  the  conventional 
methods  of  measurement  based  on  the  conditions  of  an 
unbounded  semispace  or  full  space  can  be  applied  also  in 
mine  galleries, 

(c)  to  determine  +-he  absolute  value  of  conductivity. 


3.1.1.  Inhomogeneity  of  rock 

It  could  be  shown  that  the  socalled  Wenner  arrangement  is 
well  suited  for  conductivity  measurements  in  mine  galleries. 
For  inhomogeneous  media,  measurements  following  the  Wenner 
method  yield  a  mean  value  of  conductivity.  Therefore  it  is 
possible  to  vary  the  range  covered  by  the  lines  of  force  and 
thus  also  the  range  for  which  the  mean  value  is  taken, by 
varying  the  electrode  spacing.  Measurements  at  many  different 
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siderable  electrical  inhomogeneity.  Scattering  of  the  re¬ 
sults  of  measurements  about  a  mean  value  of  conductivity  was 
the  greater,  the  smaller  the  space  covered  by  the  measurement. 
At  electrode  spacings  of  the  order  of  30  -  40  m,  however, 
good  agreement  of  all  measured  values  was  attained.  This 
means  that  the  rock  within  small  ranges  shows  considerable 
electrical  inhomogeneities,  the  mean  value  of  conductivity 
which  is  decisive  also  for  the  propagation  of  electrical 
waves  over  greater  distances,  is  highly  significant. 

3.1.2.  Influence  of  the  gallery  cavern  on  the  result  of 
measurement 

The  conventional  methods  of  determining  the  electrical 
conductivity  of  rock  require  the  measurement  to  be  made  on  a 
semi-space  or  full  space.  The  results  measured  in  a  gallery, 
however,  are  affected  by  the  gallery  itself.  The  resulting 
error  will  be  negligible  if  the  field  of  electric  flow  in  a 
gallery  agrees  with  the  undisturbed  field  within  a  large  re¬ 
gion,  Under  what  conditions  this  is  the  case  was  studied  by 
three  different  methods: 

(a)  From  an  analytical  representation  of  the  current  distri¬ 
bution  between  two  electrodes  in  a  homogeneous  full  space 
it  was  calculated  how  large  the  electrode  spacing  must 

be  so  that  only  10°/a  of  the  total  current  flows  through  a 
circle  whose  center  coincides  with  the  point  that  bisects 
the  line  connecting  the  electrodes  and  whose  plane  is  per¬ 
pendicular  to  the  connecting  line,  with  the  radius  of  the 
circle  being  equal  to  the  diameter  of  the  gallery  [Ref.  4, 
chapter  3.1.1.J  • 

(b)  The  voltage  profile  along  the  line  connecting  the  elec¬ 
trodes  and  along  a  straight  line  parallel  to  it  on  the 

opposite  side  of  the  gallery  was  measured  and  compared  with 
the  voltage  prefile  valid  for  a  5>omogeneous  full  space  and 
the  corresponding  conductivity.  The  measurement  was  made  such 
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that  a  voltage  probe  was  kept  at  M  during  the  entire  measure¬ 
ment,  the  second  one  was  connected  to  the  points  of  measure¬ 
ment  in  between. 
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Using  the  denotations  of  the  figure  above  we  obtain  for 
the  conductivity: 


(3.1) 


Prom  the  measured  voltage  values  on  the  electrode  side  the 
corresponding  conductivity  values  were  calculated  and  from 
them  the  arithmetical  mean  value  &  was  determined.  In  a  dia¬ 
gram,  the  theoretical  curve  of  the  voltage  trend  was  plotted 
for  this  mean  value,  furthermore  the  measured  voltage  trend 

v 

along  the  line  connecting  the  electrodes  and  on  the  opposite 
side  of  the  gallery  (Fig.  3*'l)*  o  is  a  suitable  scale  factor 
which  permits  a  good  comparison  of  the  measured  and  the 
calculated  results.  Prom  a  certain  electrode  spacing  onward, 
the  three  discussed  curves  agree  well  except  for  deviations 
due  to  inhomogei  eities.  If  the  value  thus  obtained  is  used 
as  a  minimum  electrode  spacing  for  conductivity  measurements, 
the  conductivity  value  is  obtained  in  good  approximation 
(  C 4  J  ,  chapter  3.1.2). 

The  statement  that  a  correct  result  for  d  is  obtained  if 
the  discussed  curves  in  the  region  restricted  by  the  probes 
are  in  good  agreement,  is  identical  with  the  statement  that 
the  ratio  of  veltage  between  the  probes  on  the  electrode  side 
and  the  voltage  on  the  opposite  side  U^/U^  must  be  unity. 
Measuring  this  ratio,  however,  we  shall  find  considerable 
deviations  from  the  value  expected  according  to  the  geometrical 
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conditions,  because  of  rock  inhomogeneities.  If  this  ratio 
is  determined  at  numerous  sites  of  measurement  for  a  certain 
probe  distance,  the  arithmetical  mean  of  all  resulting  values 
will  come  very  close  to  the  correct  value,  since  the  rock 
inhomogeneity  is  completely  irregular. 


3.1.5*  Determination  of  rock  conductivity  in  a  mine 
gallery.  Solution  of  the  respective  boundary  problem. 

The  influence  of  the  gallery  cavern  on  the  result  of  a 
measurement  by  means  of  the  Wenner  method  was  also  studied 
theoretically.  [l7j.The  results  were  in  good  agreement  with 
the  measurements . 

The  problem  to  be  solved  is  determined  by  the  following 
data: 


Ai>=  0  anywhere  in  the  medium 
§  =  const.  /  R  for  r  <=o_| 

§  %  const.  [l/R^  ~  l/RgJ  near  the  electrodes, 

=  0  on  the  cylinder  mantle  (R^,  R2  ... 
the  electrodes). 

The  resulting  solution  is  a  Fourier  series: 
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K  (ma)  =  modified  Bessel  functions, 

n\  /  i 

£,  z,vp  =  cylinder  coordinates, 

h  =  half  the  electrode  distance 

a  =  radius  of  cylinder. 

This  series  w as  computed  on  an  electronic  computer  for 
h  =  12  m  and  the  angles  0  and  f  =  t:.  The  cylinder  radius 
was  assumed  to  be  1  m.  The  calculation  yielded  the  result  that 
the  change  in  potential  trend  for  the  chosen  electrode  spacing 
is  smaller  owin g  to  the  cylinder  influence,  than  is  the  deviation 
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of  the  measured  values  from  the  calculated  potential  trend, 
due  tc  the  inhomogeneity  of  rock.  Furthermore,  the  current 
density  was  calculated  from  the  expression  for  the  potential 
and  from  it  the  portion  of  current  flowing  through  a  circle 
with  the  radius  R  in  the  plane  z  =  0  for  different  electrode 
spacings  (total  current  =  l),  A  comparison  of  the  result  with 
the  analogous  value  for  a  homogeneous  full  space  suggested 
the  conclusion  that  the  influence  of  the  cylinder  can  be 
neglected  from  an  electrode  spacing  of  2h  =  20  m  onward. 

3.1.4.  Results  of  d-c  measurements 

Measurements  of  this  kind  consistently  showed  that  the 
distance  of  probes  under  the  given  circumstances  (gallery 
diameter  =  1.5  m)  must  be  larger  than  7  -  8  m  in  order  to 
yield  a  sufficiently  accurate  result. 

Since  the  voltage  is  measured  by  means  of  a  compensating 
circuit,  the  relative  error  of  measurement  grows  smaller  the 
higher  the  voltage  to  be  measured.  Owing  to  the  great  probe 
distance  and  the  resulting  high  voltage  to  be  measured,  the 
Yfenner  method  proved  more  suitable  than  any  other  conventional 
arrangement.  The  electrical  conductivity  value  will  only  be  a 
mean  value  Decause  of  the  inhomogeneity  of  rock  throughout 
the  accessible  mines.  Here  again,  the  \7enner  method  proved 
highly  suited,  as  the  influence  of  small -dimensioned  inhomo¬ 
geneities  is  comparatively  small  because  of  the  large  probe 
distance  as  compared  to  the  electrode  spacing.  A  suitable 
choice  of  array  will  thus  yield  a  value  averaged  over  a  large 
region. 

From  the  discussed  values  of  the  U^/U^  ratio,  an  accurate 
statement  can  be  made  on  the  inhomogeneity  of  rock.  Especially 
for  small  probe  distances  the  values  of  the  Uj/Ujj  ratio 
scatter  considerably  (Fig.  3.2).  This  means  that  the  specific 
conductivity  within  a  small  region  changes  considerably,  where¬ 
as  the  average  value  determined  by  means  of  a  four-electrode 
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configuration  is  almost  constant. 

3.1.5*  Results  of  measurements  in  various  mines 

Yfith  allowance  for  the  discussed  conditions,  direct-current 
conductivity  measurements  were  made  in  numerous  mines.  The 
results  shall  be  reviewed  briefly  hereinafter. 

Bleiberg:  Shaft  Rudolf 

In  dead  rock,  values  ranging  between  0.9  and  1.1  •  10 
were  measured.  The  slight  difference  between  the  two  extreme 
values  suggests  that  the  rock  bodies  are  largely  homogeneous, 
yet  the  number  of  measurements  made  in  this  connection  is  too 
small  as  to  allow  a  final  statement  in  this  direction.  Measure¬ 
ments  in  strongly  inhomogeneous  ore  (zincblende)  showed 
greater  fluctuations  in  conductivity  than  those  in  dead  rock, 
the  values,  however,  had  the  same  order  of  magnitude.  Only 

in  comparatively  homogeneous  zincblende  on  level  10,  the  con- 

—2  —3  —1  “1 

ductivity  v/as  much  higher  (6  •  10“  to  1  *  10~^  Q  m  ). 
Salzgitter ,  mine  Konrad 

In  this  mine,  the  measurements  in  ore  as  well  as  those 

in  dead  rock  yielded  only  slight  fluctuations.  Rock  containing 

-2 

small  quantities  of  ore  yielded  values  between  3  •  10  and 
1  .  10  ^  Q  ^m  \  measurements  in  oolithic  rock  yielded 
the  value  3*1  ♦  10  ^m  \  The  maximum  deviations  from  this 
valxie  were  9%, 

Hannover,  potash  mine  of  Empelde 

Measurements  in  dry  salt  were  impossible  with  the  available 

measuring  device  because  of  the  low  conductivity  which  was 

-9  -1  -1 

smaller  than  1  »  10  J  Q  m  . 

Siegerland,  mines  of  Friedrich  Wilhelm  II  and  Georg 

The  rock  of  both  mines  is  strongly  inhomogeneous,  the 
ore  (siderite)  that  occurs  there  is  mixed  with  dead  rock. 
Measurements  in  dead  rock  and  in  ore  yielded  values  of  the  same 
order  of  magnitude?  in  both  cases,  however,  the  values  obtained 
at  different  sites  varied  considerably.  The  measured  values 
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-4 

were  uniformly  distributed  over  the  range  between  2  •  10 
-3  -1  -1 

and  8  •  10  Q  m 

Heerlen,  coal  mine  of  Oranje  -  Nassau  III' 

Only  a  small  cross  section  was  available  for  measurements 

in  coal.  Therefore  the  number  of  measurements  is  very  low. 

—8  —1  —1 

Results:  coal  in  the  direction  of  seam:  ^ =  8.0  •  10  Q  m 
coal  vertical  to  direction  of  seam: 

<o=  1.5  •  10“ W\ 

The  slate  zone  on  top  of  the  coal  from  an  electrical  view¬ 
point  is  largely  homogeneous  with  a  conductivity  of 
_i 

5  .  10  m  , 


3.2.  Frequency  dependence  of  &  and  £.  of  rock 


The  frequency  dependence  of  the  electrical  parameters 
of  rock  was  mainly  measured  on  samples,  some  measurements 
were  also  made  on  solid  rock.  The  measurements  comprised 
the  real  part  and  the  imaginary  part  of  the  complex  dielec¬ 
tric  constant: 

t*  =  £’  -  i£"  .  (3-4) 

The  real  part  is  the  relative  dielectric  constant  of 
the  medium,  whereas  the  imaginary  part  yields  the  con¬ 
ductivity  according  to  the  following  relation: 

«  £0  .  (3.5) 


3*2.1.  Measuring  device 

A  Schering-bridge  circuit  was  used  for  the  measurement. 
This  bridge  indicates  the  capacity  and  loss  factor  tan  cT  = 

=  €,"/ e  '  of  the  complex  unknown  impedance. 

For  the  measurements,  the  specimen  was  put  between  the 
plates  of  a  capacitor.  This  measuring  capacitor  can  be  re¬ 
presented  as  a  parallel  circuit  of  capacitance  and  resistance. 


Between  the  quantities  Rp  and  Cp  of  the  equivalent-circuit 
diagram  and  the  real  part  and  imaginary  part  of  the  complex 
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dielectric  constant  there  exist  the  following  relations: 

-  Cp/00 

t"  -  l/“HpC0  (5.6) 

and  US'  =  £o/RpCQ, 

-12  -1  -1 

where  £  =  8.86  •  10  a.sec*v  m  is  the  absolute  di- 

o 

electric  constant  of  the  vacuum,  is  the  capacitance  of 
the  measuring  capacitor  without  rock  samples. 

In  most  cases  it  is  more  favorable  to  conduct  the 
measurement  by  means  of  a  substitution  method,  because  of  the 
inevitable  wiring  capacitance  and  the  restricted  measuring 
range  of  the  bridge.  The  points  that  must  be  taken  into 
account  in  this  connection  and  the  technique  of  measurement 
are  described  in  detail  in  |_19l  and  foL 

3*2.2.  Review  of  possible  error  sources. 

The  quantities  of  and  show  a  great  dependence  on 
frequency  and  humidity.  Various  attempts  of  explaining  the 
frequency  dependence  have  been  described  in  publications. 
Above  all,  polarization  phenomena  on  the  electrodes  are  made 
responsible.  It  is  well-known  that  two  metallic  electrodes 
in  a  dissociate  solution  have  the  same  effect  as  a  capacitor 
of  high  capacitance,  if  a  current  is  produced  between  the 
plates.  The  most  essential  contribution  to  this  capacitance 
is  made  by  a  thin  layer  on  the  electrodes,  partly 
owing  to  a  change  in  concentration  on  the  plates  caused  by 
the  current,  and  partly  by  the  resulting  concentration 
potential.  For  direct  current,  this  capacitance  may  assume 
very  high  values  (up  to  0.02  f/cra2).  It  depends  on  the 
electrode  material  and  on  the  surface  of  the  electrodes. 

When  using  alternating  current,  these  processes  are  of 
importance  only  if  the  ions  have  time  to  concentrate  on 
the  electrodes  during  half  the  period.  In  his  theoretical 


<t  J, 
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studies,  Joliffe  [26J  found  that  the  capacitance  decreases 
with  the  root  of  the  frequency,  using  the  assumption  that  the 
concentration  potential  alone  contributes  to  the  capacitance 
due  tc  polarization.  Various  experiments  confirm  the  theo¬ 
retical  result  whereas  ethers  contradict  it.  A  satisfactory 
explanation  of  the  frequency  dependence  of  the  quantities 
and  &"  in  this  manner,  however,  is  impossible. 

In  a  report  also  dealing  with  the  measurement  of  electrical 
rock  parameters,  Tarkhov  [ 7j  states  that  the  frequency  de¬ 
pendence  of  the  quantities  f  1  and  g,"  is  not  observed  on  rock 
samples  separated  from  the  electrodes  by  thin  mica  plates, 
he  therefore  assumes  the  electrode  polarization  to  be  respon¬ 
sible  for  the  freqtiency  dependence. 

Since  the  above  report  has  been  frequently  quoted  by 
numerous  authors ,  the  Tarkhov  measurement  was  repeated  in 
order  to  get  an  explanation  of  the  influence  of  electrode 
polarization  on  the  result  of  measurement.  The  sample  was 
separated  from  the  electrodes  by  two  plastic  foils.  First, 
the  dielectric  constant  and  the  loss  factor  of  the  insulating 
material  were  determined.  For  £}  the  frequency-independent 
value  of  =  1.21  was  obtained,  the  loss  factor  being 
negligibly  small. 

For  the  capacitance  of  a  capacitor  shown  in  the  diagram, 
the  following  value  was  obtained: 


C  » 


i  »F‘ 
o 


£l’£2 


2dl62  +  d2ei 


(3-7) 


1 

i. 


f— 5 


r 
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The  fact  that  the  frequency  dependence  of  the  capacitance 
o'*  the  measuring  capacitor  is  not  observed  with  the  rock 
sample  being  insulated  from  the  electrodes  by  plastic  foils, 
led  Tarkhov  to  the  conclusion  that  the  frequency  dependence 
of  £’  in  rock  can  be  explained  by  electrode  polarization. 

This  conclusion,  however,  is  incorrect.  Under  the  assumption 
that  ££  should  be  constant,  the  expression  (3*7)  clearly 
shows  the  condition  for  which  the  frequency  dependence  of 
C  vanishes: 

2dx  8'2  »  d2  s{  .  (3.8) 

The  measurements  described  in  the  present  report  show 
that  s'  is  always  very  large  especially  in  the  range  of  a 
great  frequency  dependence  so  that  the  above  condition  (3.8) 
is  satisfied.  Thus  we  obtain 

VP* 

C  =  — 2d'~-  =  const  •  (5.9) 

Measurements  of  electrical  rock  parameters  on  samples  may 
be  affected  by  errors  due  to  possible  surface  conductivity. 
Whether  or  not  this  is  the  case  has  been  studied  by  using  a 
guard  ring  capacitor  for  the  measurement.  In  all  cases  of 
measurements  with  and  without  a  guard  ring,  the  same  results 
were  obtained;  this  means  that  there  exists  no  essential 
surfa  e  conductivity.  The  dependence  of  electrical  parameters 
on  the  frequency  thus  is  not  due  to  disturbing  effects  of 
measurement,  but  is  an  actual  property  of  rock. 

3.2.3.  Reproducibility  of  measuring  results  [l6j 

As  the  electrical  parameters  depend  on  the  frequency 
and  on  the  water  content  of  the  samples,  it  must  be  studied 
how  far  the  humidity  values  can  be  reproduced.  Since  a 
cooling  liquid  must  be  used  for  drilling  and  cutting  the 


samples,  the  water  content  frequently  varies. 

The  water  content  of  the  samples  was  determined  by 
weighing,  the  change  in  humidity  was  brought  about  by  long-time 
storage  in  an  atmosphere  of  constant  humidity.  The  entire 
reproducibility  examination  was  made  so  that  the  electrical 
parameters  were  determined  at  different  humidity  values, 
the  water  content  at  first  being  increased  stepwisely, 
beginning  with  dried  samples,  until  no  weight  increase  could 
be  observed.  Then,  the  sample  humidity  was  changed  in 
opposite  'direction,  again  several  measurements  of  the  elec¬ 
trical  quantities  being  made.  A  good  reproducibility  of  the 
electrical  parameters  was  obtained  in  dependence  of  humidity, 
but  only  in  thin  samples  (thickness  %  4  -  5  mm).  In  thicker 
samples,  the  changes  in  results  are  probably  due  to  the  fact 
that  the  homogeneity  of  samples  regarding  their  water  content 
is  not  preserved. 

The  dependence  of  the  dielectric  constant  on  the  hu¬ 
midity  decreases  as  the  frequency  increases,  approaching  a 
constant  value.  At  high  humidity  values,  the  frequency  de¬ 
pendence  is  quite  distinct,  in  dry  samples,  however,  it  is 
very  low. 

The  manner  in  which  the  conductivity  depends  on  the  fre¬ 
quency  is  hardly  affected  by  the  water  content.  The  absolute 
value  strongly  depends  on  the  humidity,  changes  by  a  factor 
of  1000  were  observed. 

3.2.4*  Explanation  of  the  frequency  dependence  by  dielectric 
polarization 

An  explanation  of  the  observed  frequency  dependence  of 
the  electrical  rock  properties  resulted  from  a  comparison 
with  K.V7.  Wagner's  theory  [28]  of  dielectric  polarization. 

The  most  simple  mathematical  model  of  dielectric  po¬ 
larization  was  given  by  Debye.  He  assumed  the  polarization 
to  consist  of  two  portions:  one  portion  which  is  assumed 


a. 
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to  have  the  same  phase  as  the  external  field  E  ,and  a  portion 
P0  with  a  time  displacement  by  the  relaxation  time  x  vnth 
respect  to  the  E  field.  This  model,  however,  is  not  sufficirt 
for  describing  the  frequency  dependence  of  electrical  data  of 
a  number  of  dielectrics  as  well  as  rock.  The  frequency  be¬ 
havior  is  represented  much  better  by  Wagner's  expansion  of 
this  theory.  For  only  one  relaxation  time  we  obtain  the 
relaxation  function 

$(t)  •  .-t/T,  (3.10) 


indicating  how  the  final  value  of  the  polarization  portion  Pg 
(not  of  like  phase  as  E)  is  attained.  Expanding  Debye's 
theory,  Wagner  assumed  that  the  behavior  of  numerous  dielectrics 
is  determined  by  an  infinite  number  of  relaxation  cimes  which 
are  distributed  about  a  dominating  relaxation  time  (z  = 

=  t/t  )  according  to  the  probability  function 

W(t)dT  =  2  dz  .  (3»ll) 


Under  this  assumption  we  obtain  for  the  relaxation  function; 

+oo 


oo 


(3.12) 


with  b  being  a  measure  of  the  width  of  the  distribution 
function  which  becomes  the  wider  the  smaller  b. 

For  the  quantities  £'  and  we  obtain; 


OO 


0 


cosh  (2b2  zq  -  l)u 
cosh  u 


du]  , 
(3.13a) 


e  - 


cosh  2b2  z  u 
_ o 

oush  u 


(3.13b) 
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with 

z  =  In  cut 
o  o 

u  =  z  +  z  =  In  u>t  . 
o 

These  equations  contain  four  constants  which  all  have  a 
certain  physical  meaning:  k  being  a  measure  of  the  change 
in  dielectric  constant  within  the  frequency  range  in  which  a 
dielectric  polarization  occurs,  6^  being  the  dielectric 
constant  of  the  medium  at  the  upper  end  of  this  frequency 
band,  b  being  a  measure  of  the  density  of  distribution  of 
relaxation  times  about  the  prevailing  value  of  tq  and  f  =  l/r0  • 
being  the  frequency  at  which  the  dielectric  losses  have  a 
maximum.  By  means  of  a  graphical  method  described  by  XI. A. 

Yager  [29]  which  has  already  been  discussed  briefly  in  [4|> 
the  results  obtained  on  rock  samples  have  been  compared  with 
.is  theory.  Furthermore,  the  dependence  of  the  constants 
contained  in  Wagner's  theory  on  the  water  content  of  the 
samples  was  studied.  A  dolomite  sample  (D  III)  from  the  mine 
of  St.  Gertraudi  whose  electrical  data  were  determined  for 
nine  different  humidity  values  was  treated  in  detail.  It  was 
four  out  that  the  studied  types  of  rock  have  very  similar 
electrical  properties,  therefore  the  results  obtained  for 
D  III  shall  be  discussed  first. 

5.2.5.  Measurements  on  D  III 

The  nine  measurements  on  the  sample  Dill  were  made  in 
two  groups:  the  measurements  1  -  6  were  made  at  low  humidity 
values.  It  has  already  been  said  that  the  water  content  was 
varied  by  storage  in  an  atmosphere  whose  absolute  humidity  was 
changed.  The  measurements  7>  8  and  9  followed  about  two 
months  later.  In  the  meantime,  the  samples  were  stored  in  a 
damp  atmosphere  (about  30$)  then  the  water  content  was  in¬ 
creased  in  the  same  manner  as  above  (measurements  7  and  8), 

In  order  to  increase  the  water  content  further,  the  sample 


OggW  ; * 


was  stored  in  tap  water  for  another  week  and  then  it  was 

superficially  dried  with  blotting  paper  (measurement  9)» 

Figs,  3*3  and  3.4  show  the  results  of  all  nine  measurements. 

Evaluation  of  the  measurements  1  -  8 . 

The  results  of  the  measurements  1-9  are  3hown  in 

table  3«1.  All  four  constants  of  Wagner's  theory  could  only  be 

evaluated  in  the  measurements  3-8,  because  in  Yager's 

report  only  a  limited  range  of  the  parameter  values  zq  has 

been  dealt  with.  An  exact  determination  of  b  and  f  was  not 

m 

possible  at  the  restricted  range  of  measurement  of  the  bridge. 

For  measurement  1,  b  was  approximately  estimated  and  for 

measurement  2  also  f  . 

m 

For  b  and  f  ,  a  unique  dependence  on  the  humidity  was 
obtained:  both  quantities  increase  with  the  vmter  content. 

This  means  that  the  frequency  at  which  the  maximum  losses  occur 
is  displaced  toward  higher  values  and  that  the  frequency 
range  in  which  losses  occur  becomes  narrower. 

Also  for  <5^  and  k  strong  changes  were  observed  owing 
to  the  change  in  the  water  content  of  the  samples,  but  a  clear 
relation  could  not  be  observed. 


Table  3.1 

Result  of  evaluation 


Nr.  of 

measurement 

weight 

Sample  D  III,  = 

d  =  6.7  am 

P  6 

=  2.65 

f 

m 

Pf 

£' 

OO 

c-k 

.  k 

OO 

exp 

1 

2 

34.9518 

34.9523 

0 

0.00143 

%  0.10 
%  0.12 

*0.017 

9 

3 

34.9531 

0.00372 

0.19 

0.36 

9.2 

256 

9.5 

27.6 

7 

34.9542 

0.00687 

0.20 

1.49 

17.8 

174 

10.2 

17.0 

8 

34.9554 

0,0103 

0.22 

9.2 

10.0 

18  p 

10.2 

17.6 

4 

34.9560 

0.0120 

0.24 

11 

11 

268 

9.6 

27.9 

5 

34-9563 

0.0129 

0.2  6 

41 

16 

235 

10 

23.5 

6 

9 

34.9569 

35.0157 

0.0146  0.29 
0.183 

81 

17 

200 

10,9 

10.4 

18.3 

p  =  percent  by  weight  of  the  water  contained  in  the  sample 
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It  is  striking  that  ^  is  larger  and  k  is  smaller  for 
the  measurements  7  and  8  than  it  is  in  the  other  measure¬ 
ments.  The  reason  for  this  fact  so  far  has  not  been  ex¬ 
plained,  yet  it  is  quite  clear  that  the  long-time  storage 
of  the  samples  used  for  the  measurements  7  and  8  caused  a 
certain  change  in  dielectric  constants,  whereas  the  constants 
which  are  decisive  for  the  frequency  trend  of  £"  remain 
unaffected  fj]  . 

The  quality  of  the  reproduction  of  experimental  data  by 
Wagner's  theory  and  the  accuracy  with  which  the  distribution 
constant  b  can  be  determined  were  discussed  in  1 4~1  •  Only 
at  very  email  frequencies,  deviations  from  theory  were 
observed:  in  general  the  values  are  in  good  agreement. 

Prom  the  good  agreement  between  the  experimental  data 
and  the  values  obtained  from  \7agner's  theory  it  may  be  con¬ 
cluded  that  the  dielectric  losses  observed  on  the  partly 
moistened  rock  sample  in  the  concerned  frequency  range  can 
be  explained  by  polarization  phenomena.  Because  of  the 
complex  nature  of  rock,  one  relaxation  time  is  not  sufficient 
for  describing  how  the  portion  of  polarization  not  of  like 
phase  as  E  is  attained,  but  for  a  satisfactory  explanation 
a  distribution  of  relaxation  times  about  a  dominating  value 
must  be  applied. 

Results  of  measurement  9 

The  result  obtained  for  measurement  9  on  samples  largely 
saturated  with  water  could  not  be  reproduced  by  Wagner's 
theory  jl9j .  Since  the  assumption  that  the  properties  of  the 
water  contained  in  the  sample  are  decisive  for  its  electrical 
properties  was  quite  evident,  and  of  water  were  also 
measured.  The  values  of  the  mixed  body  consisting  of  water 
and  rock  (measurement  9)  were  found  to  lie  between  the 
values  of  its  components.  Despite  the  small  portion  by  volume 
of  water  0.55$)  the  frequency  dependence  of  the  two 
quantities  is  similar  to  their  frequency  dependence  for  water. 
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This  fact  and  the  absence  of  a  frequency  dependence  typical 
of  Wagner’s  theory  permits  the  conclusion  that  the  occurring 
losses  are  no  longer  due  to  dielectric  polarization  if  the 
pore  volume  of  rock  is  largely  filled  with  water,  but  due  to 
line  losses,  with  the  charge  carriers  migrating  through  the 
rock  without  local  restrictions.  These  losses  are  inversely 
proportional  to  the  frequency.  In  solids  they  contribute 
essentially  to  the  total  losses  only  at  very  low  frequencies 

(  l3o]  and  [31]  ). 

3.2.6.  Measurements  on  other  rock  samples 

Studies  on  the  frequency  dependence  of  £'  and  we  re 

made  also  in  other  rock  samples:  dolomite  samples,  new  red 
sandstone,  ore-containing  rock  from  the  mines  of  Bleiberg 
in  Karnten,  Austria,  and  Fiisseberg  in  Siegerland,  Germany, 
as  well  as  samples  from  the  oolithic  ore  of  the  Konrad  mine 
at  Salzgitter. 

For  the  quantities  b  and  f  a  clear  dependence  on  the 
water  content  was  obtained  for  all  samples,  similar  to  that 
of  D  III,  not  so  for  the  quantities  and  k.  The  highest 
value  of  b  was  obtained  for  new  red  sandstone  which  has  the 
highest  porosity  of  all  examined  types  of  rock.  In  all  samples, 
an  essential  difference  regarding  the  frequency  dependence 
of  t'  and  was  observed,  depending  whether  the  sample  was 
moistened  in  parts  or  saturated  with  water.  Whereas  in  the 
first  case,  an  evaluation  of  the  measured  results  by  Wagner's 
theory  was  possible,  thus  explaining  the  principal  contribution 
to  the  dielectric  losses  by  polarization  phenomena,  this  was 
not  the  case  in  saturated  samples,  where  the  dielectric 
properties  3eem  to  be  determined  by  those  of  the  contained 
water. 

Only  the  frequency  dependence  of  the  values  measured  on 
the  sample  of  oolithic  iron  ore  could  not  be  reproduced  by 
a  relaxation  function  for  which  a  distribution  of  the  relaxation 
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times  about  a  predominating  value  of  to  is  assumed.  More¬ 
over,  the  form  of  frequency  dependence  would  suggest  the 
assumption  that  there  exist  two  different  prevailing  re¬ 
laxation  times.  An  exact  evaluation  in  this  direction,  however, 
has  not  been  attempted  so  far.  On  saturation  with  water,  the 
frequency  dependence  of  sJ  and  of  these  samples  agrees 
with  that  for  the  sample  1  III  so  that  again  the  occurring 
losses  on  the  whole  may  be  classified  as  line  losses  p4 1  • 

3.2.7.  Measurement  of  the  frequency  dependence  on  solid  rock 

On  solid  rock,  the  resistance  R  between  two  metallic 
electrodes  was  measured  as  dependent  on  the  frequency.  From 
it,  6"  was  determined  as  follows  [32I: 

•  lnf  • 

It  must  be  mentioned  that  the  contact  resistance  of  the 
electrodes  ,.ias  not  been  taken  into  account  in  the  above  ex¬ 
pression,  the  measuring  technique  thus  yielding  too  small  an 
absolute  value  of  conductivity.  The  impedance  between  the 
electrodes  was  also  measured  with  a  Schering  bridge.  However, 
only  the  real  part  was  taken  into  account,  since  the  earth 

O  , 

capacitance  of  the  measuring  device  can  hardly  be  determined, 
its  elimination  by  a  substitution  method  is  almost  impossible 
when  measuring  on  solid  rock.  For  applying  a  substitution 
method,  first  a  balanced  bridge  is  necessary  with  a  hot 
electrode  that  is  not  connected.  The  measuring  bridge  here 
lies  opposite  a  wall  with  a  uniform  potential.  Balancing 
with  a  connected  hot  electrode,  however,  leads  to  an  electric- 
potential  gradient  along  this  wall.  The  earth  capacitances 
of  the  lines  and  the  instrument  casings  therefore  have  diffe¬ 
rent  effects  on  the  result  of  measurement  in  either  case. 

Since  the  earth  capacitances  have  the  same  order  of  magnitude 
as  the  capacitance  to  be  measured,  the  latter  cannot  be  de- 


termined  with  sufiicient  accuracy.  Yet,  the  resistance  R 
can  be  determined  as  it  is  independent  of  the  capacitance 
and  of  the  alignment  of  lines. 

For  measurements  on  solid  rock,  battery-fed  instruments 
must  be  used  as  voltage  source  and  detector.  Since  one 
electrode  is  connected  with  the  bridge  screening,  the  re¬ 
sistance  of  the  electrodes  is  decreased  when  using  mains-fed 
instruments  owing  to  the  grounding  point  of  the  network 
being  connected  parallel  with  it;  thus,  the  existing  conditions 
are  not  clearly  defined. 

Results  of  measurement  on  solid  rock 

In  all  measurements,  the  dependence  of  conductivity  on 
the  frequency  was  only  small.  For  measurements  on  solid 
rock,  the  frequency  dependence  is  much  smaller  than  for 
measurements  on  saturated  samples.  Both  results,  however, 
cannot  be  reproduced  by  Vf-.gner's  theory.  At  low  frequencies, 

S’  is  consistent  with  the  direct  current  values  measured  at 
the  same  point. 

The  absolute  conductivity  value  measured  on  samples 
saturated  with  water  is  always  one  to  two  decades  lower  than 
the  values  measured  on  solid  rock.  Since  an  examination  showed 
that  the  electrical  data  of  the  water  used  for  saturation 
differ  but  slightly  from  those  of  the  water  in  the  mine, 
the  conclusion  was  drawn  from  this  deviation  that  the  original 
saturation  degree  of  the  samples  cannot  be  restored  by  sub¬ 
merging  the  samples  in  water.  A  similar  conduction  mechanism 
resulted  from  the  similarity  in  the  frequency  dependence 
of  samples  thus  saturated  and  of  solid  rock,  v/hich  differs 
from  the  mechanism  observed  for  samples  moistened  in  parts. 
T/hereas  for  the  latter  the  data  of  measurement  can  be 
well  explained  by  polarization  phenomena,  a  suitable  expla¬ 
nation  of  the  frequency  dependence  for  saturated  rock  (for 
measurements  on  samples  as  well  as  solid  rock)  can  be  given 
under  the  assumption  that  in  the  frequency  range  in  question 
line  losses  constitute  the  principal  portion  of  all  dielec- 


4.  Mathematical-physical  basis  for  the  propagation  and  ex 
citation  of  VLF  waves 


4.1.  Introduction 

It  has  already  been  mentioned  that  the  examination  of 
wave  propagation  in  an  underground  mine  is  subjected  to 
numerous  restrictions  such  as,  e.g.,  difficulties  of  trans¬ 
portation,  lack  of  space,  restricted  possibilities  of  setting 
up  antennas,  and  the  points  of  measurement,  etc.  Primarily  it 
is  of  greatest  importance  to  choose  the  suitable  type  of 
antenna  for  the  intended  studies.  The  choice  had  to  be  made 
between  linear  electric  antennas  (extended  live  wire)  or 
current  coils.  More  complex  types  of  antennas  and  antenna 
arrays  with  several  elements  were  not  suited,  partly  because 
of  the  small  financial  means  available,  and  partly  because  of 
the  above  restrictions. 

It  was  possible  to  record  only  the  near  field,  since  the 
wavelength  A.  in  the  interior  of  the  earth  at  the  used  fre¬ 
quencies  still  has  the  order  of  up  to  10^m  (depending  on  the 

conductivity),  whereas  below  ground  the  distances  of  measure- 

3 

ment  cover  only  some  10  m.  On  the  other  hand,  the  distances 
within  which  measurements  were  possible  are  large  as  compared 
to  the  dimensions  of  available  antennas  so  that  the  antenna 
at  these  sites  could  always  bo  looked  upon  as  being  a  dipole. 
For  this  reason,  our  attention  concentrated  mainly  on  the 
dipole  radiation  in  a  conducting  medium;  its  theory  was 
studied  in  detail  and  compared  to  the  experiment.  A  comparison 
of  the  performance  of  an  electric  di_ole  with  that  of  a 
magnetic  dipole  showed  the  latter  to  be  superior,  a  fact  which 
has  already  been  found  out  by  other  authors,  e.g.,  [2l]  . 

.A  theoretical  explanation  of  this  fact  is  given  in  4»2.  Let 
us  only  mention  that  the  energy  consumption  in  the  immediate 
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neighborhood  of  the  magnetic  dipole  is  smaller.  Therefore 
preferably  frame  antennas  were  used  for  exciting  a  VLF 
field  below  ground. 

In  another  section,  4*3>  the  theory  of  a  magnetic 
dipole  embedded  in  a  conducting  medium  is  described  in  de¬ 
tail.  The  well-known  Hertzian  solution  for  a  dipole  in  free 
space  is  expanded  for  a  conducting  medium  with  a  complex 
dielectric  constant  t .  Problems  regarding  the  validity  of 
the  solution  for  our  special  case,  c.g.,  allowance  for  qua- 
drupole  radiation,  embedment  of  the  dipole  in  a  hollow 
space,  etc.  are  discussed  subsequently.  Let  me,  furthermore 
mention  that  comprehensive  numerical  calculations  of  the 
magnetic  field  strength  for  wide  parameter  regions  have 
already  been  discussed  in  [l]  -  i_4_]. 

In  section  4 • 3 »  the  space  was  assumed  to  be  homogeneous 
and  unbounded.  In  4*4  the  geometry  of  the  actual  measuring 
arrangement  is  accounted  for,  and  a  review  is  given  on  the 
possibilities  of  calculating  the  radiation  field  of  mag¬ 
netic  dipoles  in  a  conducting  semi-space  having  a  ground  - 
air  interface.  The  studies  are  made  in  correspondence  with 
the  fundamental  studies  by  Sommerfeld  [22]  also  resulting 
in  a  separation  of  the  field  into  three  different  portions 
being  due  to  a  direct  wave,  a  wave  reflected  by  the  boundary 
layer,  and  a  ground  wave.  The  latter  probably  transfers 
most  of  the  energy  over  large  distances.  On  checking  the 
experiment,  an  expansion  of  the  theory  to  the  geometry  of 
a  semi-space  did  not  yield  any  noticeable  effects;  the 
space  in  the  range  of  measurement  therefore  was  assumed  to 
be  a  full  space.  Hence,  the  studies  described  in  [5]  were 
not  continued  in  1965 . 

Finally,  the  conductivity  determination  also  comprised 
a  theoretical  estimation  of  the  validity  of  the  Wenner  method 
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when  applied  to  a  geometrical  position  that  differs  from  that 
of  a  full  space  or  semi-space.  These  considerations  are 
summarized  in  4*5. 


First  we  shall  give  a  brief  reviev;  of  the  used  electro¬ 
magnetic  parameters 

iB*  =  jiu>yi(<5  +  i  w?)]  ^  a*  ficovitfj  = 

=  (^)l/2(l  +  i)  -  8(1  +  i)  .  ^  (1  +  i) 

% 

B  =  complex  propagation  constant 

6  -  conductivity 

l  =  dielectric  constant 

a)  •■=  angular  frequency 

]i  =  magnetic  permeability 

<5  =  depth  of  penetration 

A.  =  wavelength 
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H*  =  II  SL 
n*  4k 


-iB*R  - 


R 


jiB*  +  ^  sinJ^-. 


(4.2) 


For  the  magnetic  dipole: 

rr  XkX  t  I  • 

=  An  R  L1 


K1  e'iS  R  I.  11 

___  - - -  luJ£,  +  rr  + 

471  R  L  ^  iwpR2  ■ 

n 


p  _  K1  e 
"  4n 


-iB*R 


sin 


J 


R 


[iB*  +  jj  sin/$  , 


(4.5) 


(4.4) 


t 


I  =  current  intensity 
K  =  magnetic  flux 
1  =  length  of  antenna 


=  polar  coordinates 
in  space. 


The  difference  to  dipole  radiation  in  vacuo  lies  in  the 
fact  that  the  term  in  parentheses  in  Equs.(4.l)  -  (4»4) 
is  complex  and  that  furthermore  an  exponential  attenuation 
occurs.  In  contrast  to  air  where  the  antenna  may  be  looked 
upon  as  a  combination  of  its  ohmic  resistance  Rq  and  its 
radiation  resistance  R^,  we  here  have  another  energy- 
consuming  resistance,  namely  that  of  the  surrounding  medium. 
This  fact  is  of  much  greater  influence  for  an  electric 
dipole  than  it  is  for  a  magnetic  one.  Calculating  the  energy 
flux  through  a  spherical  surface  with  the  radius  R  we  obtain 
for  both  cases: 
electric  case: 


_  e  _  8 n  /.Ii  \2  e 
=  5 


-2SR 

"r 


T  .  2SB  2S  1  "1  /. 

fc*B  +  “ +  r + -3  (4,5) 


(in  vacuo  W®  =  40u2  (^)^  )> 
or  transformed  for  R  »  A  , 


W.e 


0  -4?R 

a  V/e  — - 

R»A  v  6VV 


Magnetic  case: 


4s,  p 

>2  ¥)r  [> 


A 

4kR_ 


f°r  a  vacuum). 


(4.6) 


For  R  » /\  this  yields: 


J* 

tf*  ,  =  /^e 

R»/\  v  2a 


Let  us  assume  that  there  exists  a  radius  R  so  that  the 

o 

energy  radiation  of  both  dipoles  through  a  spherical  sur¬ 
face  of  this  radius  is  the  same.  This  yields  the  following 
relation  between  the  energies  and  W®  which  thus  holds 
also  for  the  moments  K1  and  II  of  the  two  dipoles: 


yj 

V 

r 


18. 109  he^Ro^ 

f  h  (R  ) 
nr  o ' 


where  h^  and  h^  express  the  terms  in  parentheses  in 

Equs.  (4.5)  and  (4.6) ,  respectively.  Under  this  assumption, 

6 

the  relation  ¥_  for  a  second  radius  R  and  thus  approximately 
also  for  R  yield.3 : 

£  W 

an  expression  which  is  always  smaller  than  unity. 

Owing  to  ohmic  losses,  the  energy  flux  of  the  magnetic 
dipole  is  thus  weakened  by  the  surrounding  medium  less 
than  the  electric  dipole.  This  fact  is  discussed  in  detail 

in  [4]  . 


4.3.  The  field  of  the  magnetic  dipole 


The  solution  of  the  Maxwell  equations  for  a  magnetic 
dipole  at  the  origin  of  a  polar  coordinate  system  R,  Ay 


4 


4 
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yields  the  following  expression  for  the  field  components: 


n  K1  ,iB*  A  1  v  -iB~R  q 

HR  -  2?  +  e  003  ^ 

ft 


K1  / .  .1  In 

H  =  U»e.  -  Jr  +  — j) 

3.a>pR 


-iB*R 


R 


sin^-  (4.7) 


K1  /  •  1  \  © 


-ifl  R 


R/  R 


(7 

sm  <iA-  , 


Since  the  component  vanishes,  the  field  is  fcransversally 
electric  (TB) .  The  electric  lines  of  force  describe  circles 
around  the  z-axis,  whereas  the  magnetic  lines  of  force  are 
located  in  the  meridional  plane .  As  there  exists  no  ^-de¬ 
pendence*  all  meridional  planes  are  equivalent.  At  the 
beginning  of  4*2  in  connection  with  the  field  parameters, 
approximate  expressions  have  been  used  for  the  complex  pro- 
pagation  constant  0  and  for  the  complex  wave  impedance  , 
These  expressions  can  be  used  if  the  imaginary  part  prevails 
in  the  complex  dielectric  constant 


*  16" 
*  - £  ■  =&; 


so  that 


■ht 


At  a  frequency  of  5  kes  of  the  lower  boundary  of  the  VLF 
range,  this  is  the  case  from  about  0<l  5*10  '  S/m  onward. 

The  difference  between  the  expressions  (4*7)  and  the 
wave  field  in  vacuo  mainly  lies  in  the  fact  that  an  exponential 
attenuation  occurs  owing  to  the  complex  propagation  constant 


iu)6Po  /,  .  N 

13  (1  + 1)- 


Furthermore,  the  complex  terms  in  parentheses  of  (4*7)  bring 
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about  that  the  radial  energy  flux  through  a  spherical  surface 
even  in  the  immediate  neighborhood  of  the  antenna  depends 
on  the  radius  of  this  surface;  the  concept  of  radiation 
resistance  used  for  a  vacuum  thus  is  useless,  since  this 
expression  here  becomes  a  position-dependent  quantity. 

It  has  already  been  mentioned  in  4*2  that  this  is  due 
to  the  fact  that  the  energy  leaving  the  antenna  is  transferred 
toward  infinity  by  an  electromagnetic  wave  only  in  parts, 
whereas  part  of  it  is  lost  already  in  the  neighborhood  of 
the  transmitter  by  heating  of  the  medium.  The  greater  the 
distance,  the  smaller  the  loss  portion.  In  the  far  field 
it  is  given  only  by  exponential  attenuation  which,  because  of 

B  - 


isproportional  to  the  root  of  frequency  and  conductivity. 

At  a  frequency  of  3  kcs,  this  attenuation  factor  lies  between 
—4  —2  —1  *  - 

10~ +  and  10~  m  (cf,  [2|  ),  depending  on  the  quantity  of 

conductivity  (E>  =  10  ^  -  5*10  2  s/m). 

Attenuation  is  the  main  reason  also  for  a  distortion  of 
the  radiation  characteristics,  again  compared  with  the 
vacuum  case. 

The  voltage  induced  in  the  receiving  antenna  is  proportio- 
nal  to  the  absolute  magnetic  field  strength  at  the  site  of 
measurement,  therefore  the  measurable  quantity  is  the  amount 
of  the  expressions  in  (4.7)*  We  thus  obtain  the  following 
relation: 


IH 


-  (  IHjjl  2  +  |HjJ2)1/2 


„  -6R  f  . 

=  Me  |^4  cos  rj 


2  Q.f  1  2 B  2B2 

p6  +  p5  +  ,4 


+  sin2$ 


1  +2|+2Z  +  *§!  +  4§i 

R6  R5  R4  R5  R2 


1/2 


(4.8) 


M  =  moment  of  the  antenna. 


+ 


1/2 


sin 


(4.9) 


In  these  expressions,  the  displacement  current  has  already- 
been  neglected  (£,«--).  The  function  |H|,  |H  |  ,  jH^iwas 
tabularized  in  wide  parameter  ranges  as  dependent  on  the 
distance  R  and  the  angled,  being  represented  in  a  diagram. 

In  most  cases,  the  curve  trend  was  found  to  be  in  good  agree¬ 
ment  with  the  measured  curves  [ 2 J , [ 3l  • 

In  the  direct  neighborhood  of  the  antenna  where  the 
distances  of  the  measuring  arrangement  are  still  comparable 
to  the  antenna  dimensions,  the  expressions  (4.7)  are  to  be 
looked  upon  as  the  first  terns  of  an  expansion-  by  spherical 
harmonics  (cf.  f2j).  Allowing  for  another  term,  namely  the 
quadrupole  term,  we  obtain  an  additional  field  component: 


H 


=  3a  (3  cos 2  $  -  l)  e"lS  Rj— 


Q 


R 


3iB 

R5 


s*2 


r 


„Q  3sin  2 ^  -ifl*R  (  6  6iB*  3B*2  16*3  ) 

H^  =  aQ  2  *  1r4  r3  -  r2  -  R  J* 


If  the  amount  of  the  field  strength  is  calculated  with 
allowance  for  these  additional  expressions,  the  differences 
as  compared  to  pure  dipole  radiation  are  not  noticed  from  a 
distance  R  =  10  m  onward.  Hence,  quadrupole  radiation  has 
no  essential  influence  and  may  thus  be  neglected. 

The  effect  of  the  semi-space  surrounding  the  antenna 
in  the  form  of  a  secondary  field  reflected  by  the  boundary 
layer  to  the  conducting  medium  can  also  be  neglected.  In 
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accordance  with  experimental  results,  a  theoretical  examination 
(C2]ardl3j  )  shows  no  noticeable  effects  whatsoever. 

4.4.  Radiation  fields  of  magnetic  dipoles  under  the  earth's 
surface 


So  far,  mathematical  and  physical  studies  were  made 
for  a  transmitting  antenna  and  receiver  embedded  in  a  con¬ 
ducting,  homogeneous  and  unbounded  medium.  This  corresponds 
to  reality  only  so  far  as  the  distance  between  the  transmitting 
and  receiving  devices  in  most  cases  was  too  small  for  a 
reflection  by  the  air  -  soil  interface  or  the  well-known  phe¬ 
nomenon  of  a  ground  wave  being  of  any  importance.  Thus,  the 
direct  wave  was  received  whose  principal  term  is  proportional 
to  e~ ^  lS  R/r5> 

At  greater  distances  (field  strengths  were  re¬ 
corded  up  to  a  distance  of  7  km),  a  deviation  of  the  measured 
curve  from  the  trend  expected  in  accordance  with  the  theory 
of  (4.3)  was  evident.  The  possibility  of  considering  the 
geometry  of  a  homogeneous  semi-space  was  therefore  taken  into 
account  and  the  radiation  field  of  horizontal  magnetic  dipoles 
at  a  certain  depth  h  under  a  plane  surface  was  calculated  by 
applying  the  following  three  methods ; 


(a)  the  saddle  point  method  (yields  the  far  field) 

(b)  a  VLF  approximations  (yields  a  review  of  the  wave  propagation 
at  low  frequencies  and  smaller  distances) 

(c)  by  series  expansions  of  the  Sommerfeld  integrals. 

All  three  methods  are  characterized  by  the  integral  representation 
of  the  Hertzian  vector.  They  only  differ  in  the  manner  the 
integrals  were  evaluated.  According  to  Sommerfeld  [2l]  the  field 
is  reduced  to  a  two-component  vector  potential  it  =  (Ttx,0,7t2) 
by  the  relation: 

2  v 

tt  n  *-  ^  4  „  _  _  s\ 

“E  =  WE  'VE  ‘  b'La'u'  U--LV  “g 
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+  grad  div  7? 
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(the  index  E  refers  to  the  earth), 
ad  (a): 

Satisfying  the  boundary  conditions  for  z  =  0  (earth's  surface), 


we  obtain: 


iB  JR  _•  o  f 


Ex  R 


+  2  J  Ho  ^BER  sineC)  e 


iB]r(z+h)  ooeac/ 


f  («£/)  sin«6  ck. 


(4.10) 


TC„  =  cos 
Ez 


f-rik/n o1)(ber  sin^  e 


iBg(z+h)cos  <£/ 


ffU)  sin«6  dc6, 


where  f(oC)  and  g(oO)  are  certain  functions  obtained  from  the 
boundary  conditions.  The  integration  path  which  must  encircle 
a  branch  line  has  to  be  chosen  adequately  j  5]  •  A  saddle 
point  integration  f  3}  yields  the  following  expression  for 
and  for  the  component  of  the  field  strength,  with  the.  displacement 
current  being  neglected  (cf.  f  3_j  ): 
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Bc  =  constant  of  propagation  in  vacuo. 

/  \  ,  .  -iu)t  . 

The  time  factor  under  (a)  is  assumed  to  be  e  m  con¬ 
trast  to  all  other  calculations, 
ad  ( b ) 

Calculation  of  the  dipole  field  by  VLF  approximation.  The 
region  for  which  the  field  strength  expressions  derived  by 
the  saddle  point  method  are  valid,  is  somewhat  restricted, 
since  in  some  parts  of  the  ca.lculation  high  values  of  the 
dimensionless  quantity  of  BR  were  assumed,  equal  to  those 
characterizing  the  far  field  of  a  dipole.  Thus,  the  distance 
from  the  extreme  far  field  to  several  wavelengths  was  covered 
by  separate  expressions  in  order  to  on  the  one  hand  be  able 
to  calculate  the  field  strength  for  all  distances,  and  on  the 
other  hand  to  have  a  possibility  of  comparing  the  near  field 
and  the  far  field  which,  because  of  the  small  distance  covered 
by  the  measuring  transmitter,  is  of  special  interest.  Starting 
point  again  is  the  vector  potential  %  which  is  represented 
as  an  integral  over  the  eignfunction  )e-^Z: 


7t  =  M 
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V(z‘h)  +Hf  f(A)e*i‘(2+h)xo(A/j>)j  A 


OP 


nz  =  J  g(A)  e"^2!!  (Aj>)d/V  cos  y?  , 
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where  p  =  y  A/  -  Sg  • 
With  the  condition 


(— r~)2  »  i 


being  fulfilled,  B  may  be  neglected  as  compared  to  Bg  in  the 
expressions  g(A)  and  f(/L),  thus  yielding: 
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The  physical  situation  created  by  this  mathematical  approxi¬ 
mation  can  be  easily  interpreted  by  the  equations  (4.13) 
comprising  only  fields  induced  by  conduction  currents.  In 
the  atmosphere  such  currents  cannot  exist.  Neglecting  the 
displacement  currents  it  follows  that  the  region  for  which 
the  expressions  are  valid  is  restricted  to  distances  within 
which  the  conduction  currents  below  ground  were  mainly  induced 
by  primary  excitation,  not  by  an  induction  effect  of  the 
displacement  currents  in  the  atmosphere,  as  it  would  be  the 
case  at  large  BR,  i.e.«  in  the  vacuum  far  field.  The  re¬ 
sulting  expressions  are  therefore  applicable  to  the  near  field 

M  ' 

up  to  distances  of  several  vacuum  wavelengths ,  thus  supple¬ 
menting  the  expressions  of  the  saddle-point  method  in  the 
region  of  small  distances. 

The  integrals  (4.13)  can  be  reduced  to  well-known 
integral  expressions,  thus  yielding  the  following  relation 
for  the  horizontal  components  of  the  magnetic  field  strength; 


H 

y 


2_ 

3z: 

E 


+  h  W] 


M  [' 


3x3y  L  B 


-iBER  -i8ER'  -IV' 

- -  _  £ - .  +  ~  h -  [I* _ 

n  t  r  ^  .  a  \  ^  i 


B^z2  V  R' 


(4.14) 


I  (cC)  K 
o' 


4-15 


iB 


R' 


R 


E 


=  x2  +  y2  +  (z  +  h)2  oC  =  2 

/  2  2  77  77?  iBe 

®  \x  +  y  +  (2  -  h)  ^  = 


(R  -  (z  -  h)) 
2  (R  +  (z  +  h)) 


ad  (c): 

In  the  integrals  (4.12)  of  (b),  the  Hankel  functions  can 
be  substituted  for  the  Bessel  function  after  changing  the 
interval  of  integration.  Closing  the  integration  path  at 
infinity  where  a  branch  line  including  a  pole  in  the  immediate 
vicinity  of  the  branch  point  must  be  avoided,  the  integration 
can  be  reduced  to  the  following  integral  calculus: 
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Expanding  the  nominator  which  occurs  in  U,  we  finally  obtain 
the  following  expressions: 
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which  lead  to  the  field  strength  components  H  ,  H  and  H  : 

x  y  z 
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This  method  of  numerical  evaluation  whose  validity  again 
depends  on  the  neglegibility  of  displacement  currents,  is 
much  simpler  than  that  of  (b),  The  three  representations 
of  (a),  (b)  and  (c)  are  discussed  in  detail  [3]  and  are 
compared  numerically. 


In'  the  theory  of  conductivity  determination  by  the  Y/enner 
method  4 • 5  >  it  must  be  found  out  how  far  the  Wenner  arrange¬ 
ment  can  be  applied  to  the  underground  geometry,  namely  a 
cylindrical  cavity.  For  this  purpose,  the  solution  of  the 
following  boundary  problem  for  the  potential  4>  of  two 
electrodes  is  given  in  [4]  J 


4  -  35 


For  $  we  obtain  an  integral  representation  which  for  <?  =  a 
assumes  the  value  of 
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The  integrals  contained  in  this  expression  were  evaluated 
successfully  partly  by  asymptotic  representations,  partly 
by  numerical  calculation,  yielding  the  result  that  the 
expressions  which  are  based  on  the  \7enner  method  [4J  are 
valid  for  a  homogeneous  full  space  if  the  electrode  spacing 
is  sufficiently  large  as  compared  to  the  cylinder  radius  (=  mine 
gallery). 
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